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Semileptonic Decays of Heavy Flavors in QCD 



Lev A. Koyrakh, Ph.D. 
The University of Minnesota, 1996 

Supervisor: A. Vainshtein 

A model independent approach based on a generalization of the operator product expansion 
is used to describe semileptonic decays of heavy flavors. In the first part of the dissertation 
we calculate differential distributions in the inclusive semileptonic weak decays of heavy 
fiavors in QCD. In particular, the double distribution in electron energy and invariant mass 
of the lepton pair is calculated. The distributions are calculated as series in where mq 
is the heavy quark mass. All effects up to are included. 

Also calculated are the energy distribution and semileptonic decay width for the case 
of a heavy lepton in the final state. In the case of i?-meson decays, for h —f utv transitions 
the nonperturbative corrections decrease the decay rate by 6% of its perturbative value, 
while for h —f ctv they decrease it by 10%. 

Based on the results of the first part of the work, the full up to order l/rnq set of 
the OPE sum rules for the heavy fiavor transitions and radiative corrections to them up to 
order cisAqcd /''^Q a-re calculated. 

A new model is proposed for the inclusive semileptonic decays of the B mesons 
B IvXc, which is defined by the requirement that it should satisfy the QCD consistency 
conditions. These conditions are imposed in the form of the OPE sum rules. It is shown that 
under some natural assumptions the OPE sum rules provide sufficient number of constraints 
to fully determine exclusive contributions of a number of resonances in the final state of 
the decay. The proposed model can be used for experimental measurements of the heavy 
meson matrix element ji^ which describes the kinetic energy of the heavy quark inside the 
hadron, as well as for extraction of value of |Vc6| from experimental data. 
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Chapter 1 

Introduction 



The problem of theoretical understanding of semileptonic decays of heavy particles have 
long attracted theorists's attention. One of the most challenging aspects of the problem has 
to do with the strong interaction. Namely it is difficult to theoretically calculate different 
formfactors for exclusive decay channels. The problem is that the exclusive channels are 
formed by the strong interaction at large distances when the perturbation theory is not 
applicable. Phenomenological models have been employed in the situations when theoretical 
calculation from the first principles was not possible. Different models of the semileptonic 
decays of heavy quarks have been suggested [2], [28]. Those models are based on different 
sides of theoretical understanding of the problem. The ACCMM (Altarelli, Cabibbo, Corbo, 
Maiani and Martinelli [2]) model takes into account the so called primordial motion of the 
heavy quark inside the hadron and assigns Gaussian distribution of the heavy quark spatial 
momentum with the so called Fermi momentum pp being the width (dispersion) of that 
distribution. Of course, one of immediate consequences of this picture is that the lifetime 
of the meson becomes longer than in a free quark decay due to the Lorentz time dilation. 
The other, ISGW (Isgur, Scora, Grinstein, Wise [28], [47], [46]) model uses quark- antiqu ark 
potential in order to calculate the formfactors, which are obtained as overlap integrals 
between final and initial states wave functions of the so called light degrees of freedom of 
the corresponding mesons. 

In the recent years model independent predictions for the lepton energy distributions 
as well as double distributions in the lepton energy and invariant mass of the lepton pair 
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in the decays have been obtained [5], [6], [13], [38]. These predictions are based on the QCD 
operator product expansion (OPE) in which two new phenomenological parameters describ- 
ing heavy particles were introduced: matrix element of the chromomagnetic operator /^q 
which is responsible for the hyperfine mass splitting within heavy quark doublets, and fi'^ 
- average kinetic energy of the heavy quark inside the hadron. The results of the operator 
product expansion are related to measurable quantities via sum rules [49], [48], [9], [10]. One 
can also view the sum rules as model independent constraints which should be satisfied by 
any reasonable phenomenological model of the decay. In the present work an attempt is 
made to introduce a phenomenological model of semileptonic decays of B mesons which is 
consistent with the QCD in the sense that it satisfies the OPE sum rules. In the proposed 
model the final state of the decay is represented as sum of different resonance contributions 
along with continuum states. Then the resulting hadronic tensor is restrained by the OPE 
sum rules. It is shown that with limited number of resonances in the final state the con- 
straints could be solved exactly yielding predictions for different formfactors for the decay 
thus satisfying the QCD conditions. Once the constraints are solved, one can use them to 
calculate the hadronic tensor and differential distributions in the semileptonic decays. In 
this way the differential distributions in the decay could be described without ^-function 
singularities which appear in the pure OPE approach. One can also relate the parameter 
fi'^ to the quantities which could be (at least in principal) measured by experimentalists 
(such as hadronic structure functions, or at least the lepton spectrum) and in turn measure 
the yU^. This approach could be also useful in a more accurate determination of |Vc6| and 
\Vub\ from experimental data. 

In order to perform a lot of calculations done in this work a special program package 
was written in MATHEMATICA (Wolfram Research, Inc.). The package includes a small 
but effective "Lorentz calculator" intended for performing calculations and simplifications 
of expressions involving various tensor and gamma matrix structures, and a special program 
which effectively does the operator product expansion described in this work. 

The dissertation is organized as follows. In Chapter 2, based on the work [13], the 
operator product expansion is used to derive the differential distributions in semileptonic 
decays of heavy fiavors with a massless final state lepton (electron), then in Chapter 3, 
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based on the paper [33], the lepton energy distributions are derived for the case of a massive 
lepton (r-lepton) in the final state. In mostly review Chapter 4 the results of Chapter 2 
are used to derive the OPE sum rules (in that we follow works [48], [9], [10]), and the full 
list of the sum rules is calculated. In Chapter 6 a new model of semileptonic decays of 
heavy fiavors is proposed. The model is based on the QCD consistency conditions following 
from the OPE sum rules. In Chapter 6 the model is formulated in the heavy quark limit. 
Perturbative corrections to the sum rules are calculated in Chapter 5. In Chapter 7 the 
model is generalized to take into account perturbative and nonperturbative corrections to 
the OPE sum rules. Numerical analysis of the proposed model is made in Section 7.4. 
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Chapter 2 

Differential distributions in 
semileptonic decays of heavy 
flavors in QCD 

2.1 Introduction 

Differential distributions in semileptonic decays of heavy flavors are used for measurements 
of the CKM matrix elements, key phenomenological parameters of the standard model. To 
extract the CKM matrix elements from data one needs to disentangle the effects of strong 
interactions at large distances from the quark-lepton lagrangian known at short distances. 

Up to now essentially two approaches are applied to describe nonperturbative strong 
interaction effects in the inclusive weak decays: the naive parton model amended to include 
the motion of the heavy quark inside the decaying meson [2]; and the 'exclusive variant' 
based on summation of different channels, one by one [28]. Both approaches are admittedly 
model-dependent, neither their accuracy nor the connection to the fundamental parameters 
of QCD are clear a priori. Each of them needs an input from constituent quark model to 
parametrize nonperturbative effects. The latter play an especially important role in the 
form of the spectra near the endpoints. 

The need for the model-independent QCD-based predictions is apparent. Consid- 
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erable progress achieved recently in the theory of preasymptotic effects (proportional to 
powers of l/rriQ where itiq is the heavy quark mass) allows one to make these predictions. 

The theoretical construction presented in this chapter is, in a sense, a generalization 
and combination of the formalisms which are used in deep inelastic scattering and total 
cross section of e+e" annihilation. The expansion parameter in deep inelastic scattering is 
where Q is the momentum transfer. In the problem at hand the expansion parameter 
is niQ^ or, more exactly, the inverse energy released in the final hadronic state (in the rest 
frame of the decaying quark). 

In the classical problems of this type, like e+e" - annihilation, there are two alter- 
native ways to get predictions. The first approach having a solid theoretical justification in 
terms of OPE [57] is based on calculations in the euclidean domain where one can apply 
OPE. The contact with the observable quantities is made through the dispersion relations 
and in this way predictions for certain integrals are obtained. In the second approach we 
perform the calculations directly in Minkowski domain. Although formally this calcula- 
tion refers to large distances, from the first approach we know that in specific integrals 
large distance contributions drop out. Therefore the results obtained in this way, although 
not valid literally, should be understood in the sense of duality: being smeared over some 
duality interval the theoretical prediction should coincide with the smeared experimental 
curve. The inclusive weak decays will be treated within the second approach. The averaging 
mainly refers to the invariant mass of the inclusive hadronic state produced in the decay 
considered. 

If the invariant mass of the final hadronic state is large this is not a constraint at 
all since the theory 'itself takes care of the averaging required by duality. In the opposite 
limit, near a spectral endpoint, the smearing is not provided for free. The boundary of 
the distribution corresponds to low momentum of the quark produced (low momentum 
of the hadronic final state). At this point the OPE blows up, therefore we do not have 
any specific prediction for the distributions near the boundary. Nevertheless the integrals 
taken over the domain from the kinematical boundary up to a new boundary, defined by 
the requirement that the OPE is convergent, are predicted. In particular this integration 
domain should include the resonances range (when itiq is large a parametrically stronger 
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limitation is imposed by the fact that the heavy quark and meson masses are different). 
An example of the safe integration is the total decay width where the integration domain 
is maximal. 

Although we explicitly work in the Minkowski kinematics we always keep in mind 
the relationship to the euclidean domain and the corresponding operator product expan- 
sion. The first analysis of this type has been outlined in [54] for inclusive heavy fiavor decay 
rates. A general analysis of the semileptonic inclusive spectra along this line is presented in 
Ref. [18]. In that work it was observed, in particular, that the leading operator and those 
appearing at the next-to-leading order have a gap in dimensions of two units, and, conse- 
quently, the O^niQ^) term should be absent in certain quantities. The analysis presented 
in [18] was not backed up, however, by concrete calculations of the preasymptotic effects. 
Recently this formalism has been systematically developed and applied to the non-leptonic 
decays of heavy fiavors [6, 15] and the charged-lepton energy spectrum in the semileptonic 
decays [7] (see also [5]). 

We generalize the results of Ref. [7] to find the complete inclusive distributions in the 
semileptonic decays. The leptonic variables — Ee,q^ and qo, where Eg is the charged lepton 
energy and q is the momentum of the lepton pair ^ — are kept fixed which automatically 
fixes the invariant mass of the inclusive hadronic state. Integrating over go we obtain the 
double spectral distribution in Eg and q^. 

At the first stage we construct the transition operator T(Q X ^ Q) describing 
the forward scattering amplitude of the heavy quark Q on a weak current. Our focus is 
the infiuence of the 'soft' modes (background fields) on the transition operator T^j, which is 
expressed as an infinite series in the local operators built from gluon and quark fields and 
bilinear in Q, Q. 

The local operators are ordered according to their dimensions; the coefficient func- 
tions contain the corresponding powers of l/mg (or 1/Eh, where Eh is the energy released 
into the hadronic system). At sufficiently large nig or Eh the operators with the lowest 
dimensions dominate, and the infinite series can be truncated. Generically, we will refer 
to the power expansion as l/mg expansion, although strictly speaking it is an expansion 
^The charged massless lepton produced will be generically called 'electron' hereafter. 
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in 1/ Eh- At the next stage the matrix elements of the relevant operators over the initial 
heavy hadron Hq must be evaluated. Unfortunately, in the present-day QCD the matrix 
elements over the hadronic states are not theoretically calculable. In some instances they 
can be related, through heavy quark symmetries, to measurable quantities [24], [28]; in other 
cases they have to be parametrized. These parameters play the role analogous to the gluon 
condensate [49]. As a matter of fact, at the level of the leading preasymptotic corrections 
only two operators are relevant. The matrix element of the first one can be related to the 
mass splittings of the vector and pseudoscalar heavy mesons. The matrix element of the 
second one has the meaning of the average square of the spatial momentum of the heavy 
quark Q in Hq and the state must be treated as a parameter. 

Finally, the observed decay rates and spectra are obtained by taking the discontinuity 
of the hadronic tensor i^H q\T ^j^\H q) and convoluting the result with the lepton currents and 
appropriate kinematic factors. 

In this work we consider the differential distributions in the semileptonic decays at 
the level of 0{m'^). The differential distributions are measured experimentally in the B 
meson decays and will be used for more precise determination of Vub, for example. This 
was a primary motivation for our investigation. We would like to make it as close to the 
fundamental QCD as possible. 

The organization of this chapter is as follows. In Section 2.2 we describe the kinemat- 
ics and in Section 2.3 we present the operator product expansion. In Section 2.4 we derive 
the differential distributions for the massless lepton case. Section 2.5 is devoted to analysis 
of our distributions and limitations on the their use. Then we apply our results to deriva- 
tion of the heavy lepton energy distribution and semileptonic width in decays Hf, tvX. 
Our results are summarized in Section 3.3. Appendix B contains expressions for hadronic 
invariant functions. 
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2.2 Kinematical analysis 

We will consider the inclusive weak decays of the mesons (or baryons) with the open heavy 
flavor into the lepton pair plus (inclusive) hadronic state 

Hq{ph) Kpi) + ^iPv) + hadrons. 
Our goal in this chapter is to calculate the differential decay rate 

2.f 



dEedq'^dqo ' 

where Eg is the energy of the emitted electron and = -\- is the 4-momentum of the 
lepton pair. In order to flnd the differential distributions we need to know the amplitude of 
the process, which is given by the expression 

M = V,Q-^eV,v{X\j^HQ). (2.2) 

Here VqQ is the corresponding Cabibbo-Kobayashi-Maskawa matrix element, = qT is 
the electroweak currents, T ^ = 7^(1 + 75). (Although our theory is general we will keep 
in mind the h c and h u decays, so that Q = b and q = c or u). The differential 
distributions we are interested in are given by the modulus squared of the amplitude (2.2) 
summed over the flnal hadronic states. 

The modulus squared of the amplitude summed over the flnal hadronic states can 
be written as 

\M\^ = \VgQ\^GlMH^F^W^,, (2.3) 
where Mhq is the mass of hadron Hq , l^^j, is the hadronic tensor 

W,. = J2(27r)U\pH^ _ ^_p^)_i_(i7Q|4(0)|X)(X|j;(0)|i7Q), (2.4) 

and Z^'' is the lepton tensor 

l>^- = 8[ {Perip^r + iPenp.Y - a^^iPe ■ p.) + ^e''''"''(per(p.f]. (2.5) 

Let us introduce the hadronic structure functions Wi and parametrize the hadronic 
tensor in the following way: 

Wf,^ = -vji g^i, + vj2Vi_,v^ -i VJ3 e^^^afS VaQfS + Qf, Qu + ws (qu Vf, + q^ v^). (2.6) 



Here = (pe + Pj/)^ is the 4 - momentum of the lepton pair, = (phq)ii/Mhq is the 
4- velocity of the initial hadron (not that of the Q-quark). Note that we have omitted the 
structure g^fj/ — ^j/f^ which can not appear because of the T-invariance. The structure 
functions vji depend on two invariant variables, q ■ v and q^. In the rest frame of Hq which 
will be used throughout the paper q ■ v = qo, so Wi = Wi(qo,q^). The convolution of l^^j, 
with the lepton tensor (2.5) is given by the expression: 

W^, F'' = 4{2q^wr + [4Ee{qo- Ee)-q^]w2 + 2q^2Ee- go) W3 }. (2.7) 

We see that only three structure functions are relevant for the processes we are considering 
in this paper. At this step we encounter the third variable, the electron energy Eg = 
Pe ' PHq/Mhq, entering through the leptonic tensor. 

Finally the formulae for the differential width takes the form 

This expression concludes the kinematical analysis. Our task is, of course, the calculation 
of the structure functions Wi(qo,q^). We will proceed to this calculation in the next section. 

The phase space of the decay is described in details in Appendix A. Here let us 
write down expressions for different Lorentz invariant phase space integrations which we 
will need in the future. 



LIPS = III dEJq^dqo 



2Mg Mg-2Ee 2Mg 

dEe j dq^ j dqo 



{Mg-Mof 2Mi V 

,2 



dq / dqo / dE, 







10 -\/ In 



j dq^ j dqo j dE^. (2.9) 



gp-lql 

2 
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The limits of integration are easily established by considering kinematical boundaries 
for the particles taking part in the decay. 

We now can perform the integration over the electron energy and obtain distribu- 
tion in lepton pair energy and momentum (the distribution is in fact the integrand in the 
following expression for the total width): 



^ = \v^q\'t^J ^q'/^ diomd-e)yH + \y^2}. (2.10) 

Id tt^ Jo J\fa^ 
Integrating equation (2.8) over go we obtain the double distribution in the electron 
energy and the invariant mass of the lepton pair: 



(lE^dq^ ' 32 TT 
where 



.max 



7 = l^gQp^jTT^ / rfgo{2g'wi + [4^e(go-^e)-g']«^2 + 2g'(2^e-go)«^3}, 

^ < } TT^ ' mm 
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= {Ml-Ml + q^)l2MB. (2.12) 
Let us notice that the last distribution (2.11) contains weighted integrals of Wj's over 

go: 

^ ^max 

qQW,{qo,q^)dqo, 



in other words, moments of the structure functions, where g™™ and q^^^ are kinematical 
boundaries of go. This is important, because it is the moments of the structure functions 
which appear in the inclusive distributions are the quantities which could be reliably pre- 
dicted in QCD. We will show this in Section 4.1 devoted to derivation of the OPE sum 
rules. 

As it is well known, not all values of go correspond to physical processes. The ones 
which do lie on the physical cuts of the forward scattering amplitude T^j, (which will be 
introduced in the next section) in the complex plain go. The picture of the cuts is shown 
on Fig. 2.1. This picture could be readily understood. First, all the cuts lie on the real 
axis for only real go are physical. Different cuts correspond to different physical processes. 



B meson decay processes correspond to < go < Mb — \J Mf^ -\- q^, this cut will be called 
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Figure 2.1: Cuts of the forward scattering amplitude in the complex plane go- 

cuti. There is also a cut corresponding to collision of B and D mesons, which lies in 
Mb + \J + q2 < go < +00, this cut will be important in considering operator product 
expansion sum rules, we call it cut2. All negative go are physical, they correspond to a 1^ 
boson hitting B meson with formation of some final hadronic state X, we call it cut^. Let 
us note that separation between the decay cuti and collision cut2 is 2^ + >> Aqcd- 

2.3 Operator product expansion 

In this section we will discuss the derivation of the tensor Wi_tiy. The Operator Product 
Expansion (OPE) is similar to that in the deep inelastic scattering. It is convenient to 
introduce the hadronic tensor /i^j, (forward scattering amplitude) as follows: 

K. = ^J d'xe-^^^j^{HQ\T{j+{x)jM}\HQ). (2.13) 

Q 

The absorptive part of this tensor reduces to l^^j, discussed above 

PF^, = (l/i)discV. (2.14) 

Here disc(/i^j,) is the discontinuity of the forward scattering amplitude /i^j, on the physical 
cut in the complex plane of the variable go- Of course, /i^j, can be expanded into the same 
set of structures as l^^j, (see. Eq. (2.6)) 

h^qi^qu + hsiq^Vf, + qi_,v^), (2.15) 
11 



Figure 2.2: The tree diagram determining the transition operator T^j, in the leading approx- 
imation. The dashed lines correspond to the weak currents, the solid internal line describes 
the propagation of the quark q and the bold external lines represent the heavy quark Q. 

and the relation (2.14) implies that 

w, = 2lmh,. (2.16) 
Let us remind that /i^j, is the matrix element of the transition operator T^j, 
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V = WT^iHQ\T,AHQ), (2.17) 



T,, = i / d\e-'^-T{]l{x)],{Q)], (2.18) 



so below we will construct OPE for the product of currents in equation (2.18). Having in 
mind the relationship to euclidean analysis discussed above we will treat our expansion in 
the same way as a normal euclidean OPE. In the asymptotic limit uiq oo the hadronic 
tensor /i^j, is given by the tree graph of Fig. 2.2. This graph defines the matrix element of 
the transition operator T^j, over the heavy quark state, 

{Q\T,AQ) = -u-qT, ^_ ]_^^ r. ^Q- (2.19) 

The latter expression represents nothing else but the free quark decay. In the asymp- 
totic regime uiq oo the interaction of the heavy quark with the gluon/light quark 
medium, as well as its intrinsic motion inside the hadron can be neglected. Then 

P^ = Po^ = mQV^, (2.20) 
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where f ^ is the 4- velocity of hadron Hq . 

Equation (2.19) allows one to immediately write down the operator form in the 
approximation at hand (only operators bilinear m Q, Q are considered, see the discussion 
of other operators at the end of this section): 

- 1 2 . - 

T^i/ = — Qr^-T Fjy Q = ~7T2 ^[s'a/i^i/ + 9ai/k^ — g^i/ka — i^fn/afskp] Q"/"(l + 75)^7 

(2.21) 

where k = Pq — q. 

As we see, the two operators QjaQ and QjalsQ showed up in the operator expansion 
at the level considered. Note, that the QjalsQ term vanishes after averaging over the 
unpolarized hadronic states. 

In this paper the perturbative corrections in are not touched upon at all. As for 
nonperturbative corrections they appear due to interactions with the soft medium of the 
light cloud in Hq. By taking these interactions into account we isolate two types of effects. 
First, the fast quark q produced does not propagate as a free one, but interacts with the 
background fields; these corrections will be included explicitly into the OPE coefficients. 
Second, the heavy quark Q also does not live in the empty space; it is surrounded by the 
light cloud. In particular, due to this fact the heavy quark momentum does not coincide 
with rnqv^. This large distance effect will not be calculated explicitly, but implicitly it will 
be refiected in the Hq matrix elements of the operators in T^j,. This is in full analogy with 
what people usually do in deep inelastic scattering. The infiuence of the background field 
on the transition operator is summarized by the following expression 



T^. = - J dxe-'^^Q{x)T^S,{x, 0) E, g(0), (2.22) 

where Sq(x, 0) is the propagator of the quark q in an external gluon field A^. It is convenient 
to use the Schwinger technique of treating the motion in an external field (for a review of 
QCD adaptation see , e.g. Ref. [43]). Within that formalism the propagator Sq is presented 
by the following expression 

5',(^,0) = (a;|^^|0). (2.23) 

Here ^ = 7^(p^ + A^(X)), = g A^^T^ is the gluon field in the matrix representation. 
Furthermore, the operator of coordinate and momentum are introduced, (thus the 
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field Ai_i(X) becomes an operator function of X^), with the commutation relations 

[p^,X,] = ig^,, [X^,X,] = 0, [p^,p,] = 0. (2.24) 

The states \x) are the eigenstates of the operator X^, = 
Combining equations (2.22) and (2.23), we arrive at 

T^. = - [ dxe-'^-(x\Q(x)T^—^ Q(X)\0). (2.25) 

J jl-' - niq 

As we have discussed above the operator contains a large mechanical part (-Po)^* = 
rnqv^] the deviation from Pq will be separated explicitly 

= (Po)^ + vr^ (2.26) 

and we will expand in vr^. In this paper we will limit ourselves to the terms up to C'(vr^) 
corresponding to I/ttiq corrections. The master formulae to perform the expansion is 

r^. = - / MMQ(x)T,j—j^^^^-^ g(x)|o). (2.27) 

There is a subtle point in the description of the formalism given above. Technically 
in the computation the A^(a;) is assumed to be a c-number background field while in the 
final expression for local operators it should be understood as a second quantized operator. 
Since we are not considering any loop corrections this substitution is justified. 

Let us now discuss the set of the operators relevant to the order 0{m'^). Without 
loss of generality we can work in the rest frame of the hadron i/g, i.e. = (1,0,0,0). Only 
those operators will be retained which produce non- vanishing results after being averaged 
over Hq. The leading operator, as it was discussed above, is 

QToQ, (2.28) 

its matrix element is fixed by the vector current conservation, 

^ 'Hq\Q^oQ\Hq) = 1. (2.29) 



Equation (2.29) is given in relativistic normalization we are using throughout this paper. 
In the non-relativistic normalization there is no need in the factor 1/2Mhq in the LHS. 
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As it has been noted in Ref. [18] there are no operators of dimension 4 in the problem 
at hand. The set includes two operators of dimension 5: 

Og = '^Q<T"f'GapQ, (2.30) 

0^ = -QT>^Q = Q7v^Q, (2.31) 

where a"f^ = ^ (j"^//^ — 7'^7"), and Gap = 9 G'^^^jT"' is the gluon field strength tensor. The 
classification above takes into account the fact that the quark field Q satisfies the equation 
of motion. In particular, it stems that the operator QQ is not independent but is reducible 
to three operators (2.28), (2.30) and (2.31): 

QQ = QioQ-7r-TQ^'Q + 

Q'j"'^Gaf3Q + O(ma^). (2.32) 



4 ml 

To get Eq. (2.32) we observe that the lower component of Q is related to the upper one in 
the following way 

and the difference between QQ and QjoQ is due to the product of the lower components. 
(Here and below we will stick to the Hq rest frame.) 

A few other useful relations which can be obtained in the same manner and are valid 
at the level 0{m'^) are: 

QlT^Q = ^Q{t^^ -'-(jG)Q + 0{mQ^), (2.34) 

QlT^loQ = 0{m-Q^), (2.35) 
gvrog = ^g(7r2 - ^(TG')g + C(mQ2). (2.36) 

A few comments are in order here concerning the actual technique of constructing the 
OPE. Since we work in the Hq rest frame it is convenient to compute different components of 
T^v separately. Too, Toi, Tjo and Tij. The calculation itself is a straightforward although 
rather tedious procedure of expanding the denominator in Eq. (2.27) in /vr using the 
properties of the 7 matrices, the commutation relation 

= iG^i, (2.37) 
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and equations (2.34) - (2.36). 

Notice that we must keep the terms of the first order in ttq and of the second order 
in TT, since 

TToQ =^-^^Q + 0(mQ'). (2.38) 
Next, observe that the Green function in the background field can be written as follows: 
1 . ^ , ^ 1 



1 

n 



(Z'-^ + m,)-. (2.39) 



To transpose l/II with Fj, it is convenient to use the identity 

ir„ = r„i + i[r„.nii 

= r4 + i[r,„i.Gli (2.40) 

Acting on Q and using the equations of motion we can now substitute l/II in both terms 
on the right-hand side by 

^ (2.41) 



uIq — rn^ — 2Vq + ' 

-2^ 



provided that we limit ourselves to terms up to 0(mQ ). The second term in (2.40) can be 
simplified even further since here we can additionally neglect vr in P = Pq + vr. 

We split the calculation into three parts: Vector X Vector, AxialxAxial and Vectorx 
Axial in correspondence with the structure of T ^ as a sum of vector and axial vector, 
= 7^ + 7^75. The full hadronic tensor /i^j, is given then by the following expression: 

^t^i^ = + + ^i^iT + h"^^^ = JJJ^i^Q l^/T/ + T/^if + + rj/l-ffg) • (2.42) 

The complete expressions for the hadronic invariant functions are given in the Appendix B. 
In the order 0{mQ^) they are defined by the matrix elements of operators Ofj, O-,^ given 
by eqs.(2.30), (2.31): 

^^(^qIQ \ <r,.G''Q\HQ) = ^±-{B\b{a ■ B)h\) = (2.43) 
where B is the chromo-magnetic field operator, and 

^ ■Hq\Q7v'Q\Hq) = i,1. (2.44) 
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The parameter /Jq coincides with m^jj introduced in [14]. For mesonic states it is expressible 
in terms of the quantity measured experimentally - the hyperfine mass splittings, and it 
has the zero value for baryonic states of the type of Aq; the parameter has the meaning 
of the average square of spatial momentum of the heavy quark Q in the hadronic state Hq. 
The two quantities /^q and often appear in the combination ~ I^g^ equation (2.32). 

The value of /^q is known from the hyperfine mass splitting between the members 
of the heavy quark symmetry doublet B and B* , for it one has the following expression 

= ^(M|. - M|) ^ 0.36 GeV^. (2.45) 

As for the value of it is not known that well. Some QCD sum rules estimates 
of were made in [4] Also the inequality > /^q was proven in [52]. Let us say here, 
that fi'^ is an important phenomenological parameter which enters different characteristics 
of the B mesons and should be measured in experiments. One of the goals of this paper is 
to formulate constraints on the phenomenological structure functions which follow from the 
OPE sum rules and then fit them using the measured lepton energy spectrum. In this way 
we hope to be able to extract the QCD parameter fi'^ from the experimentally measured 
lepton energy spectrum. Its value is expected to be around 0.5GeV^. 

Let us now introduce some definitions for the quark masses. In the works The heavy 
hadron mass is related to the heavy quark mass in the following way: 

MH,=mQ+A + 4^ + 0{^, ^), (2.46) 

where Q is a heavy quark ( in this paper it is 5 or c), Hq is the lowest lying meson containing 
Q, A is the so-called "mass" of the light degrees of freedom . The last expression is in fact 
an expansion of the hadron mass in l/mq. The fi'^/2mQ part of the hadron mass is due 
to the motion of the heavy quark inside the hadron, fiQ/2mQ is due to chromomagnetic 
interaction between the heavy quark and light degrees of freedom, A is the itiq independent 
part of the heavy meson mass. 

Let us also introduce the following notation: 

^5 = MB-m5 = A+ ^' and = Mz, - = A + ~ f ^ (2.47) 
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Figure 2.3: The tree diagram determining the operator without the heavy quark Q. Now 
the bold internal line describes the propagation of the heavy quark Q and the solid external 
lines represent the quark q. 

then 

_ = _ + _ ,^,,), _i. _ _ M + o, _L, i^), ,2.49) 

and 



^^ = iA-^)i^-^nOi.^.^, (2.50) 



The last comment of this section is about the operators which are not bilinear in Q, 
Q fields. The simplest example of appearance of such operators is given by the diagram of 
Fig. 2. 3 where the heavy quark Q propagates between the current vertices. This diagram 
is similar to the one of Fig. 2.2, and the corresponding operator follows from Eq.(2.21) by 
substitution Q =^ g, niq =^ mg, q^. The additional term in T^j, has the form 



i ^ 2_ 

fi- niQ [q^ - 

(2.511 



The matrix element of the operator qj"q over the Hq state counts the number 
of quarks q and is not small in general. The operator coefficient given by Eq.(2.51) is 
particularly large when q^ rriQ. 

In terms of the intermediate hadronic states in the forward scattering off Hq this 
contribution is due to states in crossing channel containing two Q quarks - the problem was 
pointed out in Ref. [18]. The cross-channel is not related to the weak inclusive decays under 
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consideration. It is reasonable to accept the duality between the operators without heavy 
quark fields and the cross-channel contributions having nothing to do with heavy fiavor 
decays. We rely on the assumption that we can consistently omit the crossing channel 
together with operators in T^j, related to this channel. 

2.4 Calculation of the differential distributions for massless 
lepton in the final state 

The differential distributions we are interested in are determined by equation (2.8) con- 
taining three structure functions wi, W2 and w^. They are obtained from the results for 
hi (see Appendix B) by taking the imaginary parts of the corresponding functions (see. 
Eq. (2.16)). The imaginary parts are due to the poles of hi and are obtained through the 
relations: 

1 (-l)*^ d"^-^ , , , 

Im— = 7r-^ -Siz), 2.52 

where z is given by 

z = niQ - 2mQ qo + q"^ - m^. (2.53) 

We do not present here the expression for the triple differential distribution which can be 
easily obtained by combining the equations (2.8), (2.16), (2.52) and expressions for hi from 
the Appendix B. 

Although the result is derived for the physical quantity d'^T / dEedq^dqo, it can not 
be directly compared with the experimental data. An obvious signal for this is the presence 
of the delta function and its derivatives. It is not surprising because we are sitting now 
right on mass-shell of the g-quark. As we discussed in the introduction our results should 
be understood in the sense of duality: that is that the predictions should be smeared 
over certain duality interval. At the moment we have no purely theoretical tools to fix 
the size of the duality interval, therefore we are forced to rely on qualitative arguments 
and experimental data. For example the duality interval for go can be inferred from the 
distribution in the invariant mass of the final hadronic states. Our ^-functions refiect the 
resonance structure at low invariant masses. The smearing interval should be chosen in 
such a way as to cover the entire resonance domain up to the onset of the smooth behavior. 
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Instead of smearing of the distribution one can calculate the average characteristics like the 
total width T or (M^), where Mx is the invariant mass of the final hadronic states. The 
power corrections we have calculated will enter in a specific way in each particular quantity. 

Now let us proceed to the calculation of the double differential distribution d^T/dEg dq^. 
To this end we must integrate over qo, rather simple exercise with ^-functions. However if 
one would perform the integration by merely substituting 

"mn + q'^ — rnl 

10 - = (2.54) 

2 niQ 

and taking the derivatives in the case of 6' and 6" , one would get the wrong answer. The 
point is that integration domain in go has a boundary from below 

go > i^e + (2.55) 

which corresponds to AEgE,^ > q^. Therefore one should take into account the fact that 
go can not cross the boundary (2.55). For that we introduce d(qo — Eg — q^ /AE^) into the 
integrand. The occurrence of the ^-function is important for the integration of 8'{qo — ^q) 
and 8"{qo — ^o) which leads to appearance of ^(^q — Eg — q^ /AEg) and 6'{qQ — Eg — q^ /AEg) 
in the double distribution d^T / dq^ dEg because of differentiation of the ^-function . The 
final formulae for the double differential distribution in the lepton energy Eg and q^ takes 
the form: 



dxdt ' "'^l 96 7r3 



\^^q\ x'{6il- t)il -P-X + xt)+ 

Gq [l-5p + 2t + lOpt + lOxt - lOxt^ - 



{-1 + 6p - 5p'^ + X - 5px + t - 2pt + Sp'^t + xt+ ISpxt + 
5xh - 2xf - lOpxt^ - lOx^t^ + 5x^t^) S{{1 - t){l -x)- p)] + 
Kq [-3 + ?,p + At- Apt - 6xt + Axf - 
{I - 2p + p^ - "ix + 3px -3t + 2pt + p^t + llxt - 3pxt - 
3xH - Qxt^ - 2pxt^ + 2xH^ + xH^) - t){l - x) - p) + 
{I - p - x + xt){l - t){l -2p + 

p^ - 2xt - 2pxt + xh^) S'{{1 - t){l - x)- p)]}. (2.56) 
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Here we have introduced the dimensionless variables 



x = 2EelmQ, t = q^/2mQEe, (2.57) 

and the parameters 

p = my ml, Gq = Mc/m^, Kq = ^l/m^. (2.58) 

Let us emphasize that the scale mg used in equation (2.57) is the heavy quark mass and 
does not coincide with Mhq which is normally used in the experimental distributions. 

The fact that OPE generates corrections only of the order of 0{mQ^) (terms pro- 
portional to Kq and Gq) is valid for the distributions only if we use mg as a scale, i.e. in 
the variables x,t. Of course one can easily rescale them to Mhq', then the corrections of 
the order of 0{mQ^) will show up for trivial kinematical reasons. 

We can proceed further and obtain the energy spectrum by integrating over q^. The 
range of integration is given by 

0<t<l —. (2.59) 

1 — X 

The result for the energy spectrum coincides with that obtained in [7]. For the sake of 
completeness we present it here ^: 
dV \V,Q\^Glml 



dx 1927r= 



;i _ ^ _ p)2x\(l - f Yil + 2/)(2 -x) + (l- f Y(l - x) + 



5 If) 

(1 - /)[(! - /)(2 + -a; - 2/ + -fx) - i-i2x + /(12 - 12a; + 5x'))]Gq 



^(1 - f )Hl + 2f )x + ^(1 - f )(10x - 8x') + 4(3 - 4/)(2a;2 - x^) ] Aq}(2.60) 

6 p p 

where / = p/(l — x). Finally, performing the last integration over x in the domain 

<x <1- p, (2.61) 
we arrive to the total width coinciding with that in [5]: 

r = \V,q\'^^^^[^o (l + l (Gq - Kq)) -2zrGQl (2.62) 
where = 1 — 8/> + 8/>"^ — p"^ — 12/>^ log p and zi = {1 — p)"^. 



^Let us draw the reader's attention to the difference of notation: y in [7] is equai to 
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Now let us discuss the characteristic features of the double distribution (2.56). The 
most striking one is the presence of the singular terms. The technical reason for occurrence 
of those terms was that we expanded the denominator of the pole expression (2.27) in vr 
and a G. Physically this expansion reflects the shifts of the masses of particles due to the 
nonperturbative effects. As it was mentioned above these singularities reflect the structure 
of the resonance domain and the predictions suitable for comparison with the experimental 
data require smearing over the corresponding domain. To illustrate the most salient features 
of our prediction let us concentrate on the physically interesting case of the h u transition. 

For massless u quark the kinematical region of h quark semileptonic decay is shown 
on Fig. 2.4. It has the form of a square with the side equal to I in the plane {x = 
2Ee/mb, t = l2mfjEe). The right-hand side of the square corresponds to the maximal 
energy of electron Eg = mb/2 while the upper side is a maximal energy of neutrino. In the 
real B meson decay the kinematical region is certainly wider; if one neglects the pion mass 
the region is the square with the side x^ax = tmax = MB/rrif,. The origin of this window 
is related to the motion of the heavy b quark inside the B meson. In our calculations we 
account for nonzero momentum of the b quark in the form of expansion which produced 
singular S and S' terms on the boundary. It is possible to show (see refs.[7],[3I]) that the 
expansion breaks down at distances ~ (Mb — mb)/mb near the boundary, so we need to 
integrate our distributions over a range of the order of the window between quark and 
hadron boundaries. It is interesting to note that the distribution spreads off the distances 
of the order (Mb — mb)/mb while the corrections to integrals are only of the second order 
in l/rrif,. 

Another effect we need to account for is the structure of the resonance region near 
the low end of the hadronic invariant masses. To imitate the effect let us imagine that this 
region corresponds to the u quark fragmentation into the hadronic states with s (the square 
of the invariant mass) from s = 0tos = so = 2 GeV^. The curve corresponding to s = sq 
on Fig. 2.4 is given by the equation: 



and the resonance region should be included as a whole into the process of integration; we 
can predict the integral but not the structure. 
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2.5 Application to the analysis of the experimental data 

Our theoretical prediction (2.56) depends on the following parameters: VqQ, rriQ, nig, Kq, 
Gq. Let us remind that in this paper we do not consider perturbative in corrections (see 
Ref. [17]), which, of course, should be added. The Cabibbo-Kobayashi-Maskawa matrix 
element VqQ does not effect the form of the differential distribution; the total semileptonic 
width is proportional to |VgQp. The quark masses enter at the level of leading approxima- 
tion while Kq and Gq determine I/ttiq corrections. It is important that our differential 
distributions by themselves could be used to fit these parameters. In particular it is a good 
place to extract the heavy quark mass. 

Our purpose here is to give an idea of how important the I/ttiq corrections are in 
the case of charmless i?-meson decays {h u transition). To this end we will use the 
approximate values for the parameters m^, Kf, and Gh obtained from other sources. First, 
we use rrifj ~ 4.8 GeV as deduced from the QCD sum rules analysis of the Ypsilon system 
[51], and m„ = 0. Parameter Gh can be extracted from the B,B* mass splitting [14]: 

Gb = ^{M^{B*) - M^{B))/ml ~ 0.017. (2.64) 

As a representative value we use for the parameter Kf, the value ~ 0.02. A close value was 
obtained in Ref. [4] from the QCD sum rules. Earlier QCD sum rule result [39] was a factor 
of two higher. Notice that the sensitivity of our results to the value of Kf, is essentially less 
than that to Gb- For example, Eq. (2.62) for b u transition contains Gb + \Kb- 

In accordance with the discussion at the end of the previous section the comparison 
with experiment should include integration of our distribution (2.56) over the domain which 
includes the area adjacent to the kinematical boundary. We will choose this area to be given 
by the resonance domain (see Eq. (2.63)) with sq = 2 GeV^. 

Let us introduce the quantity: 

P(x,,Q = ^ff dxdt-f^ (2.65) 

where Xc, tc is the point in (x, t) plane sitting not too close to the boundary (outside 
the resonance range), Fq = \Vub\'^ G'^m^ / and the area of integration A{xc,tc) shown 

on Fig. 2.4 as shaded includes the resonance domain plus domain x > Xc,t > tc- For the 
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Figure 2.4: The kinematical region of the decay for h u decays in coordinates x = 2Ee/mb 
and t = l2mfjEe. The solid lines are the kinematical boundary for the h quark decay 
{xmax = 't-max = 1) ^ud the dashed lines are the boundary for B meson decay {xmax = 
tmax = MbIttHj). The area of integration for the distribution P(x,t) is shaded. It includes 
integration over the resonance domain. 

experimental distribution the range of integration should be extended to include the window 
between quark and hadron kinematical boundaries. Notice that in the limit of large itiq 
the size of the window (Mhq — mQ )/mQ is parametrically larger then the resonance range 
so/itiq. In the case of 5-quark they are numerically close. 

The function P(x,t) is plotted as a function oft on Fig. 2.5 for three values of a; equal 
to 0.3, 0.6, 0.8. The last value of x is close to the border of the resonance region beyond 
which we cannot make reliable predictions for the distributions considered. The dashed 
lines on Fig. 2.5 describe the leading order distributions in t while the solid lines include 
QCD corrections we have calculated. As we can see it from the curves, the corrections are 
negative and their relative magnitude is larger near the endpoints of the spectra. 
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Figure 2.5: The integrated distribution P(x,t) for the case b ^ u is plotted as function of 
t = l2mfjEe for few values of a; = 2Ee/mb. The dashed lines correspond to the leading 
order distribution while the solid lines account for nonperturbative corrections (see Eq. 
(2.56)). The lines stop at the border of resonance region. It follows from the picture that 
the corrections are negative. 
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Chapter 3 

Corrections to the heavy lepton 
energy distribution in the inclusive 

decays Hf^ ruX 

In this part of our work we investigate the power corrections to the heavy lepton energy 
distribution in the inclusive decays 

Hb TvX, 

where Hf, is a hadron containing heavy quark h and r is the r-lepton. We include non- 
perturbative corrections up to the order l/ml, and show that they cause the decay rate 
of i?-mesons to decrease by 6 to 10 percent of its perturbative value depending on the 
mass of the quark in the final state. The lepton mass does not effect neither hadronic nor 
weak leptonic tensors. Therefore we can use expressions (2.2)-(2.8) and (B.7)-(B.ll) for 
the hadronic invariant functions to get the corresponding matrix element. What differs the 
decay with a heavy lepton from the decay with a massless lepton is the phase space of the 
particles in the final state. 

The first experimental observations of the decay were made recently using the missing 
energy tag [1],[44]. The missing energy was associated with two z/^ in the decay chain 
b t~VtX, t~ VtX' , which made it difficult to reconstruct. The branching ratio 
was found to be 4.08 ± 0.76 ± 0.62% in [I] and 2.76 ± 0.47 ± 0.43% in [44] (more recent 
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analysis), which is compatible with the Standard Model. However because of the difficulties 
in identification of the decay mode Hf, rz/X, the accuracy of the measurements is still 
insufficient for direct comparison of the results of this paper with the data. 



3.1 Heavy lepton energy distributions 

As was discussed earlier in the semileptonic inclusive decays there are three independent 
kinematical variables: go, and g^, where Ej- is the energy of the emitted lepton (r in this 
paper), go is the energy and is the invariant mass of the lepton pair. The full differential 
distribution can be written as follows: 

^ \M{E,,q^,dq^)\\ (3.f) 



dE^ dqo dq^ f 28 vr^ Mh^ 

The matrix element of the process Ai{ET,qQ,dq^) involves the same hadronic and leptonic 
tensors as used in [f3] and which are given by the expressions (2.2)-(2.8) and (B.7)-(B.f f ). 
The corresponding phase space and the boundaries for the kinematical variables for the 
massive lepton case are briefiy discussed in Appendix A. 

The double differential distribution in Ej- and q^ could be obtained from (3.f) by 
integrating over go- The corresponding expressions are somewhat cumbersome because of 
the complicated limits of integration over go in (g^, q^) plane. Also there is little use of the 
double distribution because of the difficulties in observing the considered decays at present 
time. Therefore we will not write out the formulae for the double distributions in (g^,_E^), 
but will rather integrate Eq. (3.f) twice and get the energy distribution for the charged 
lepton in the final state. One more remark is in place here. Although we are looking at a 
hadron decay, in our dynamical approach we are considering a heavy quark decaying in the 
external field created by its interactions with the light degrees of freedom. That is why in 
the case at hand we have to use the quark kinematical boundaries rather then the hadronic 
ones. This means that in the formulas (A.7)-(A.fO) we have to use and niq instead of 
Mhi, and Mo correspondingly. We will use such a kinematical boundary in the following 
calculations. Our results then should be understood in the sense of duality (see discussion 
in Chapter 2, [f3]). 
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To get the energy distribution we integrate (3.1) over dq^ and dqo. In the dimen- 
sionless variables: 

X = , p = — , and pr = (3.2) 

nib rn^ rn^ 

the result of the integration then takes the following form: 

— — = 2 ^ X'^ - A p^ {i (I - Pq + p^) X - 2 x"^ - A p^ + Pq(?> - Pq + ?> Pt-) f + 
{-2- Pq+Qp,- PqP,- 2 p^,)f +2{l- p,f f + 

, bx"^ lApr 2 , , , , 2 2 ,o 

J^b[ ^ + - 3 (3 - + 3 />0 / + (1 + 3 /Og - 4 + - + />2) /2 + 

-{Q + Pq-Qp^+Apqp^-Qpl+Pqpl + Qp^^)— + 

3 (1 - (-1 - 2/>g + 2/>. - 2/>,/>, - />2) Z! + 4(1 - 4] + 



5 14 5 

Gfc [ 2 a; + - + 4 - y + ( -2 + 3 - - - 2 + 5 />0 / + 
-2 - ^ + 4 + 8 - ^ />2 - 2 _ 10 ^ 

6 6 Pq 



5 25 f-^ 

(-1 + />0 (3 + :^ /Og - 8 + - + 5 />2) ^ + 5 (1 - — ] }, (3.3) 

6 6 Pq Pq 



where 



and 



Pq 



(3.4) 



Glml 



The quantities Kf, and Gh are the hadronic matrix elements introduced in Eq. (2.44). 

The energy distribution (3.3) spans in x from x = 2 ^f])^ to a; = \ ^ p^ — pq. It 
includes the nonperturbative corrections - terms proportional to li'^ and Gh- The part of 
Eq.(3.3) without corrections coincides with the electron spectrum in p decay with a massive 
from Ref. [50]. In the limit />g ^ we encounter the familiar end-point singularities of 
the lepton spectrum. To get the limit right we make the following substitutions: 

Pq n-1 

and 



n>l, (3.6) 



r S'(l + Pr-X) , , 

Pq (ra-l)(ra-2) 
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and only then take pq to zero. To see that this procedure gives the right limit, it is sufficient 
to compare the results of integration of the both sides of (3.6) and (3.7) in the limits 
2 ^J])^ 1^ X < 1 + Pt — Pq, multiplied by an arbitrary integrable function and check that the 
both sides are equal to each other. The distribution for pq ^ takes the form: 
1 (IT, 



To dx 



= ^a;2 -4^^{(6 + 6/>^)a;-4a;2 - 8 p^- 



^ + ni + ..-)(-- + ^-2.^ 



3 3 '^3 3 ' 3 ' 
Gb[ f^ + S{l + p,-x){-- + 9pr-7pl + -^) + 4x + -^]}. (3.8) 

Now we can integrate distribution (3.3) to get the decay width including 
corrections: 

r = \Vq,\' ^^{(1 + ^^Y^ ) MP,.Pr) - 2 G, zr{pq,p.) }, (3.9) 

where: 

ZoiPq^Pr) = VX{1 - 7 Pq - 7 P^ + P^q - 7 Pt - 7 pI + pI + PqPri'i-'^ - 7 Pq - 7 Pt)) + 

12 pI (1 - pI) log + 12 pI (1 - pI) log |i±^ (3.10) 

Zi(pq,pr) = VX((l-pqf + (l-prf-l-pqPri4-7pq-7pr)) + 
A = X{l,Pq, p^) = 1 + p^g + pI - 2pq - 2pr - 2pqPr, 

Vq and Vt are the maximal velocities of the quark and the r-lepton produced in the decay: 

Vx 

= 5 = . • (3.12) 

^ + Pq- Pt I- Pq+ Pt 

There exists a simple relation between the two functions zo(pq,pT) and ziipq^p^^) [7] : 

/ ^ o dZo(Pq,Pr) „ dZo(pq,Pr) , . , . 

Zl{Pq,Pr) = -2pq 2 p^ +4Zo{Pq,Pr) ■ (3.13) 

Width (3.9) is symmetrical function of pq and Pt and therefore its limit when pq ^ 

T,^^o = \Vq,f^^{{l + ^^^^){l-8pr + 8pl-pt-l2pl logp,)-2G,{l-prr} 

(3.14) 

is the same as limit when Pt ^ with substitution pq =^ 

29 



3.2 Numerical estimates and experimental predictions 

To make numerical estimates in this section we will consider the i?-meson decays B tvX . 
We choose the following values for the parameters entering expressions (3.3), (3.8), (3.9) 
and (3.14) (see discussion in [7] and in chapter 2, [13]): = 4.8 GeV, rric = 1.4 GeV, 
m„ = 0, rriT = 1.78 GeV, Gh = 0.02 and Kf, = 0.017 (note that for A5 we have Cfc = and 
the corrections are much smaller). 

For the width of the decay we see that the nonperturbative corrections are negative. 
For h UTV transitions they decrease the width by 6% of its perturbative value, while 
for the h ctv case they decrease it by 10% (note that for the massless lepton in the 
final state these numbers are 4% and 5% correspondingly). The quantity insensitive to the 
uncertainties of known values of Vqh is T{h Ti'X)/T(b evX): 



The first factor of last equation describes the phase space suppression while the second 
contains the nonperturbative corrections. The corrections reduce Tc by 4% and r„ by 2% 
of their perturbative values. Note that in the quantity Tc the relative contribution of the 
corrections is almost independent of the uncertainties in the value of c-quark mass. 

The obtained energy distributions (3.3) and (3.8) could be applied to the decays 
involving the transitions h c and h u. The analysis of applicability of the distributions 
was made in Chapter 2 ([7] and [13]). There it was shown that the proper quantity to 
confront with experiment is 



This quantity does not contain the end-point singularities and is suitable for direct com- 
parison with the experimental data for x not too close to its maximal value (so that the 
operator product expansion is still valid and we are not in the resonance region). As we 
mentioned above, the hadronic kinematical region is different from the quark one. There- 
fore to compare our results with experiment the range of integration of the experimen- 
tal distribution should include the window between the quark and hadronic boundaries 
1 -\- Pr — Pq '^k X < MbIttHj -\- p^ — pg. Correspondingly, the reliable prediction can only 



T(b TvX) _ zoipg.pr) 
T(b^epX) ~ zo{pq,0) 



[l-2Gb{ 



Zl{Pq,Pr) _ ZljPq.Q) 



(3.15) 




(3.16) 



30 



0.2 0.4 0.6 0.8 1 

X 



Figure 3.1: The energy spectrum of r is plotted for b ctv transitions. The solid line 
shows the distribution with the nonperturbative corrections, while the dashed line - without 
them. For comparison, on the same plot we show the electron energy distribution for h cev 
transitions {p^ = 0). The graph can only be trusted for x < x^ax ~ 0.95. 

be made for 2y^^ < a; < x^ax, where x^ax = ^ Pt — Pq — (Mb — 'mb)/'mb. For M-quark 
Xmax ~ 1-05, for c-quark x^ax ~ 0.95. 

The lepton energy spectrum is plotted on the Fig. 3.1 for b ctv and on the 
Fig. 3.2 for h utv. The delta- functions of Eq.(3.8) are not shown on the graph. For 
comparison, on the same plots we show the energy distributions for electrons in 5 ^ cev 
and h uev transitions correspondingly = 0). 

The function 7(2;) is plotted on Fig. 3.3 for the case h ctv and on Fig. 3.4 for 
h UTV. The solid line shows 7(2;) with the nonperturbative corrections while the dashed 
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Figure 3.2: The energy spectrum (3.8) of r is plotted for b utv transitions {pq = 0). The 
solid line shows the distribution with the nonperturbative corrections, while the dashed line 
- without them. For comparison, on the same plot we show the electron energy distribution 
for h uev transitions = 0). The graph can only be trusted for x < x^^x ~ 1.05. 
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Figure 3.3: The function 7(2;) plotted for the case b ctv. The solid line shows 7(2;) with 
the nonperturbative corrections while the dashed line - without them. The graph can only 
be trusted for x < x^^x ~ 0.95. 
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Figure 3.4: The function 7(2;) for the case b utv (pg = 0). The solid line shows 7(2;) 
with the nonperturbative corrections while the dashed line - without them. The graph can 
only be trusted for x < x^ax ~ 1.05. 
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line - without them. 

Although the functions are plotted for the whole range of a;, 2y^^ < x < p^ — pq, 
we can only trust the graphs for x < x^ax- 

3.3 Summary of results of chapters 2 and 3 

Let us now summarize our results. Model independent approach to nonperturbative effects 
(f/mg)" is used for calculations of differential distributions. The effects are most pro- 
nounced near the endpoints of the spectra. We discussed how the comparison with exper- 
iment should be formulated accounting for the boundary effects. Somewhat disappointing 
is that we cannot use our results to improve an extraction of Vub by the consideration of 
dependence. Indeed, experimentally the signal of 5 ^ m is due to the range of electron 
energy Eg near the upper end where 5 ^ c is absent. However as it follows from Fig. 2. 4 
the distribution in at such energies is concentrated in the resonance domain, and no 
model-independent prediction emerges. 

We also investigated the nonperturbative corrections up to order f/m^ to the decays 
Hfj TvX and found their contribution to the r-lepton energy spectrum and the width of 
the decay. In the case of i?-meson decays the corrections could be up to 10% of the width. 
Unfortunately at present time the experimental measurements are not accurate enough in 
order to compare our results with the data. 
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Chapter 4 



Sum rules for heavy flavor 
transitions 

4.1 Cuts and sums. OPE sum rules and predictability of 
moments and inclusive distributions 

One of the goals of this work is to formulate a model for the semileptonic decays of heavy 
flavors which could describe transitions into individual hadronic final states. We will see 
that QCD provides sufficient number of constraints on the formfactors of such transitions. 
Those constraints are imposed in the form of sum rules which follow from the OPE analysis 
of the inclusive decays developed in the previous chapters. We consider these sum rules to 
be the QCD conditions on the structure functions of different final state hadronic resonances 
which should be satisfied by any phenomenological model. This is in fact what we will call 
the QCD constrained model of the semileptonic decay of the B mesons. So, we define our 
model by the requirement that it satisfies these sum rules. 

The sum rules for semileptonic heavy fiavor transitions were considered in the works 
[9], [48], [10], [11], [35], [52]. These sum rules are based on the dispersion relations which follow 
from analyticity and unitarity of the scattering matrix. This approach follows many ideas 
of the classical works [49]. Despite the fact that we are not going to use Borel transform in 
our analysis in this work, it is clear that it also could be utilized and we are going to do it 
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somewhere. 



In this chapter we are going to briefly review and summarize the results for the OPE 
sum rules derived in the above cited papers and provide the full list of the sum rules for 
heavy flavor transitions up to the order I/ttiq. This is the highest order we can go based 
on the structure functions calculations done in Chapter 2 and listed in the Appendix B (see 
also [13] and [38]). 

To be speciflc, we will be talking about h c transitions, however no sum rules will 
contain the c-quark mass in the denominator, thus they could be used for h u transitions 
as well. As we will see, it is only Ec (or _E„, correspondingly) that should remain large 
parameter for the sum rules to be valid. 

The hadronic invariant functions introduced in the previous section satisfy the fol- 
lowing dispersion relations 



Here cut means the physical cut in the go complex plain, go is some point away from the 
cut. We will distinguish two kinds of cuts: hadronic for which we keep notation cut and 
quark, which we will denote as qcut. 

In Chapter 2 (see also [13] and [38]) the tensor /i^j, was calculated by means of the 
operator product expansion (OPE) for the correlator of currents in (2.13). Let us call /i^j, 
calculated in this way hf^°'', which means that these functions vave been calculated using 
theoretical metods. The expressions for the functions obtained by means of the OPE 

are given in the Appendix B (see also [13]). For them the dispersion relation (4.1) is only 
valid for the values of go for which the OPE is valid: far from the cuts. 

On the other hand, we can represent the hadronic tensor W^,^ by inserting the 
complete set of physical states in the Eq. (2.4). We will call the tensor calculated in such 
a way W^^"^, phenomenological. In the points where the OPE is valid, the /^^^^(gcq^) 
calculated using (4.1) for wj^™ and the h^^J'°^[qo, q^) should be the same with the accuracy 
with which the OPE was calculated. Then in a remote from the cuts point go one obtains 
the following equation: 




(4.1) 




phen 



(4.2) 
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At the same time, we can also write: 

/,f-(go, q') = - ' . ^ dqo. (4.3) 

27r Jqcut qo - qo 

Now equating RHS's of (4.2) and (4.3) we have: 

1 f «)f'"(go,q')^. 1 f <n?~o,q')^. 

— / — dqo ~ — / — ago- (4-4) 

27r Jcut qo - qo 27r Jqcut % " ^0 

The point go in which the comparison of the theoretical and phenomenological in- 
variant functions is performed will be called reference pount. 

Let us describe the analytical structure of the hadronic tensor /i^j,(go,q^) in the 
complex plane of go at fixed q^. This tensor has discontinuities on the physical cuts , which 
represent possible physical processes. The semileptonic decay cut spans from go = to 



qQ = Mb — y Mjj + q^, we will call it cuti. Another cut corresponding to fusion of B and 



D mesons, goes from go = Mb + y Af|) + q^ to go = oo, we will call it cut2. There is the 
third cut going from go = — oo to go = 0, which corresponds to the process of D absorbing 
a T^-boson and going into B. This cuts structure is shown on Fig. 4.1. 




CUt\ 

decay cut 



m„-Jm ^+cf 

« V com ^ 



Figure 4.1: Cuts of the forward scattering amplitude in the complex plane go and position 
of the resonance region and continuum spectrum. 

Let us introduce a new variable 

e = MB- ^JmI, + q2 - go, (4.5) 

which describes the distance from the right end of the decay cut to the reference point go. 
Then e = Mb — M'^t + — go is the final state excitation energy counted from the Ed*. 
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Note, that in the definition of e we used the hadronic masses, and not the quark ones. Then 
for hadronic invariant functions we have the following expression: 



q') ^ - / ^ / ^ d~e + - / -^^A^dL (4.6) 

^TT Jcuti € — e Zir Jcut2+cuts € — € 



In the Eq. (4.6) e is in fact excitation energy of the final hadronic state counted 
from the D* energy. The equation (4.6) in fact is the starting point for the derivation of the 
sum rules. [49], see [10], [11], [9], [48], [52], [35]. Despite of the different ways used for deriving 
the sum rules in the original papers, they all could be obtained in a systematic way from 
the equation (4.6) (see [10]). 

The choice of the reference point e is based on the following consideration. The only 
experimentally observed processes lie on the decay part of the physical cut of the forward 
scattering amplitude. Therefore, the integral in the RHS of the Eq. (4.2) is known only for 
the part of the cut corresponding to the B meson decays. This dictates the choice of the 
reference point close to the decay part of the cut (which is cuti) and as far from the other 
cuts as possible. At the same time e can not be too close to cuti for the OPE to be valid. 
Therefore we choose the reference point in between cuti and cut2 but as close to cuti as 
possible. This will lead to the fact that in Eq. (4.6) only integral over cuti will be relevant 
for derivation of the sum rules. Note, that the value of e at the reference point is negative. 

One of important assumtions we are going to make is the assumption of local duality 
starting at the scale of M^ual — M^*. In practical terms this duality means that hadronic 
tensor Wj^™(go,q^) is the same as the theoretical one w*^''°^(go, q^)- 

The choice of the reference point close to cuti enables us to leave only integrations 
over cuti and neglect contributions of other cuts. Then with the help of the local duality 
we can write: 



-de ^ — / , de (4.7) 



where Cdual = Mg — y AfJ^^j + q^. The unusually looking lower limit of integration in the 
RHS is due to the fact that integration in quark go goes from to nif, — Ec, and, therefore, 
in e it starts from (Mb — mf,) — (Ed* — Ec). 
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For e lying on cuti where |e| ^ |e|, e < let us expand the denominators of the 
integrals over the cuti in the following series: 

= - = y-Tzr^ 4.8 

For the other cuts, on the contrary, we have different condition satisfied: |e| ^ |e|, and the 
expansion goes in the positive powers of e and not in negative. 

Now by equating the coefficients in front of we arrive at the OPE sum rules: 

/ wf^^i,, q2)e d,= ^ wf^^^e, e)e'de. (4.9) 

It is clear now that the operator product expansion approach provides answers for 
the moments of physical structure functions, but not for the functions themselves. This is 
sufficient for calculation of inclusive quantities such as inclusive differential distributions 
and total widths. In the next chapters however, it will be shown that if one consideres the 
moments (or the corresponding sum rules) as constraints on a model of the decay, then 
determination of exclusive quantities also becomes possible. 

Let us notice that the first three sum rules corresponding to k = 0,1,2 could be 
calculated using the results for /ij's derived in the work [13] (see Appendix B), since in 
that work /ij's were calculated only to the order I/itiq. The last term on the RHS in fact 
contains gluon corrections to the forward scattering amplitude and therefore all corrections 
of order due to virtual and real gluons should be included. The virtual gluon corrections 
have been calculated in the previous works [40]. Corrections due to emission of a real gluon 
will be calculated later in this work. 

4.2 Complete list of the OPE sum rules up to order l/niQ 

The RHS of the sum rules (4.9) for A; = 0, 1, 2 for j = 1, 2, 3 were calculated in the work [10]. 
Here we reproduce the results of that work and add the results of calculations for j = 4, 5. 
As it was explained in the previous section, for the moments 

rJiJ2,k 1 Phen/ 2x > 

Ij =7r % (e,q )e de 
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we have theoretical predictions. Here jj^'-''^'^ is the moment corresponding to the meson 
average of r{/i/2}- As it was explained in the previous section, for these quantities we have 
theoretical predictions. In the derivation of the hadronic structure functions the expressions 
for AA current could be obtained from the expressions for VV current by substitution 
nic —Trie. For the moments of the structure functions however this is not true because 
of the shift of variable go to e (4.5). The practical way of calculation was expansion of the 
corresponding hadronic invariant functions in 1/e (see [10] for details). In this derivation 
we used mass relations (2.47-2.50). Here we provide here the complete list of moments of 
all structure functions without perturbative corrections. The perturbative corrections to 
the sum rules will be calculated in the next chapter. 
For the VV currents: 



A. (4.10) 



4 E^ rrib^ 6 E^^ nib 4 E^ 




rVV,l ^ _^ ^ /i "^c W 1 , ( 1 \ I^G 

1 ~ o \ r 



2 V EJ AEc\ El \ EJ \ Ej / Qmb 



' ■ (4.11) 



12 Ec Qmb 4EJJ E, 



^vv,o _ , f -5 m nic \ 2 / 5 _ 2 nib _ nib nic^ \ ^ 

2 ~ E,^\QE,mb^ 2Ej'^ 2Ej')^''^UE,mb 3 EJ" 2i?/i^"' 



j-VV,l 

-'2 




(4.14) 
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jVV,2 J. 2 nib ( ( mc^\ 2 / . -, 



/r'^ = /r^ = :7^ + :7T^, (4.17) 



1 Ec\ jifj I -I Ec m/\ fii 



2 QmbJ Ey \ 4 Qmb iE^"^ Ef 



For the AA currents: 



jAA,o _ 1 ( -nic mc nic^ \ ^ 

^ ~ 2 2 I 12 4 ^^52 + 6 E^^ nib / 



-rrir nir'^ nir^ 



111 (4.22) 



4 Ec mb^ 6 Ei;^ mb 4 



(4.20) 



(l2£: + ^^ + Ji7)<l-lt>''i' •".23) 

jAA,2 ^ ^2 c c — c _^ _ _^ — E <L^^ „2 24) 

^ 1 2 2Ec 2Ec^ 2Ec^ \Q QEc 6^2 q Ec^ T 
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/2^^'° = ^+(:— + r^-Tr7^) ^+1:— (4.25) 




jVA,i _ (-^ mc\ ^ 



(4.26) 
(4.27) 



l3^'' = - t:^ - I ^ + I t^, (4.28) 



2^, V Ej 4E^^\ E^) 6E,mb 
Ec , rric^ \ III 




6 nib 4 Ec^ / E^ 



2 ' 



(4.29) 



All moments that are not written are equal to zero. Let us notice that all second 
moments are related to each other by the heavy quark symmetry, this happens because 
these moments are only known in the leading order and nonperturbative corrections to 
them have not been calculated. 

For the V — A current we have 

j-V-A,k _ jVV,k _ ^jVA,k , jAA,k 
^3 ~ ^3 ^^3 + ^3 

For the reader's convenience let us reproduce them here: 

/|'-^'° = l + ^^^^MzZfi), (4.31) 



'iE-'mb 




(4.32) 
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Ec J 3Ec 



3 ' 3^ ' E/ / 



y_A,2 T22"^6 f-. 2mb m^^ \ ^ 



SEcrrib \ 6 3m6 2Ej / 



rV-A,2 



(4.35) 



(4.39) 



/I'-^'O = 2/1^^'°, (4.41) 

j]--^'! = 2/^'', (4.42) 

ij-^^'' = 0, (4.43) 

ij-^'^ = 21^'', (4.44) 

= 21^''. (4.46) 

The listed moments could be viewed as expanded in the inverse heavy quark mass, they con- 
tain leading order contributions plus nonperturbative corrections (except second moments 
for which only the leading order is calculated). Radiative corrections to the sum rules will 
be calculated in Chapter 5. 
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Chapter 5 

Perturbative corrections to the 

sum rules 



5.1 Introduction 

In the discussion of the operator product expansion we did not take into account any 
perturbative corrections to the forward scattering amplitude. This also means that we have 
neglected contributions of perturbative corrections in the theoretical part of the sum rules. 

There are two places the perturbative corrections are coming from: virtual gluon 
exchange and real gluon emission. Therefore one can write for the structure functions: 

where Wj°^(e, q^) is the leading order in result including nonperturbative corrections 
obtained in Chapter 2, w™*(e, q^) and vj^^^i^, q^) are perturbative corrections coming from 
virtual and real gluons correspondingly. Now expresstions for each moment of each structure 
function will contain two additional parts, due to virtual and real gluons. 

The contribution of virtual gluons to the hadronic tensor could be found in the work 
[40], for readers convenience and to make this text self contained, we review the results of 
that work in the following Section 5.2. The real gluon contributions to the sum rules are 
calculated in Section 5.3. 
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5.2 Virtual gluons contributions 

In this section the results of Ref. [40] of calculations of the virtual gluon contributions to 
the hadronic vertex are used to calculate the corresponding contributions to the hadronic 
tensor and the OPE sum rules. 

The virtual gluon exchange results in the modification of the interaction vertex [40] 
(see also the review paper [41]): 

^Qb"D = (l + «'')7^ + ^''^" + cV^ (5.1) 

Tq^o = (1 + 0^)7^7^ + h'^vi'-i^ + c^^;'^7^ (5.2) 

where 

a' = —[In — In— + -{J{w)±r{w) + g{z,w)% (5.3) 

6^'^ = -'^-^[2r{w)^l + h\'^{z,w)l (5.4) 

c^'^ = T^/^r'^(^,-), (5.5) 

where 2; = mc/mb, w = v ■ v' = Ec/nic = \/l + q^/m^, and f' is the velocity of the final 
c-quark, A is the gluon mass which is needed to regularize the infrared divergency. This 
divergency is to be cancelled in the sum rules by the corresponding divergency in the gluon 
emission part of perturbative corrections. The upper signs refer to the vector current, 
whereas the lower signs refer to the axial current. Expressions for the functions /(w), r(w), 
g(z^w)^ h^'^(z^w)^ and h^2^\z^w) are given in the works [40], [41], we cite them here for 
the reader's convenience. 

w± = If ± y^w"^ — 1, (5-6) 

r(w) = I - ln(w \/ — 1), (5-8) 

f(w) = wr(w) - 2 - ^_ [^2(1 - wl) + (w'' - iy(w)], (5.9) 

Vtf^ — 1 

g(w) = ^==[L2(1 - zwl) - L2(l - zwl)] - 



\J 'Up' — \ 



z 



1 — 1'WZ + z"^ 
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[(w^ - l)r(w) ^{w-z) In z\, (5.10) 



, ln(l -t) , , 

dt^ ^, (5.11) 

t 

hl'^{z,w) = _ ^ ^ ^^^^ {2{w T 1)^(1 ±z)\nz- 
{{w ± 1) — 2w{2w 

^[In^-li^l, (5.12) 

hX^'^iz.w) = /i^'^(^-\w)- 2r(w)± 1. (5.13) 

After plugging this vertex into the expressions for the hadronic tensor, one can 
read off the corrections to the corresponding hadronic structure functions. The corrected 
expressions look as follows (nonperturbative corrections are not included): 



vjX^ = 27r(l + 2a^)(mb - - qo)S{ml - 2mbqo + - - ml), (5.14) 

= 47r(m5 + 2a^ nib + b^(mi, + nic) + 

c^(— + ■mb))S(ml - 2mbqo + qo - - ml), (5.15) 
mc 

= 0, (5.16) 



wj^ = 4tt — S(ml - 2mbqo + q^ - q^ - ml), (5.17) 
m^ 

= -27r((l + 5^ + c^) + 2(a^ + c^^))^(m2 - 2mbqo + ql - - m^), (5.18) 

mQ 

and for the V A currents only hX,^ is not equal to zero, 

= 27r(l + + a'^)8{ml - 2mbqo + q^ - q^ - ml). (5.19) 

For the AA currents expressions are obtained from the ones for VV currents by replacing 
mc —mc- 

In the last expressions we can rearrange the argument of the delta function, since in 
the sum rules we integrate over e = Mb — E^* — qo and only over the decay cut: 

S{ml - 2mbqo + qo - q^ - ^l) TT^K'^nb -E^- qo) = TT^K'^ - (Mb - mb) + (Ed* - E^)), 



where Ei = ^Mf-\-q^. The corrected sum rules are now obtained by multiplying Eq. 
(5.14)-(5.19) by the appropriate power of e, dividing by 27r and by 2Ec, dropping the delta 
functions and substituting go mb — Ec- Let us note that these corrections contain the 
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infrared logarithmic divergencies which are to be cancelled by the corresponding real gluon 
contributions. Let us call perturbative corrections to the moments including only virtual 
gluon contributions Sj. Then for the corrections to the moments we get: 

Sr^' = ^{E,-mMMB-m,)-{En*-E,)f, (5.20) 
Ec 

cVV,k ^ V I i.V/ , \ , 

^2 = 7rL2a rrib + b (mb + mc) + 
Ec 

c^(^ + m5)]((MB - nib) - (Ed* - E,))\ (5.21) 
Sl^'^ = 0, (5.22) 
Sr'' = ^{{MB-mb)-{ED*-E,))\ (5.23) 
SV'' = -:^[(6^ + c^) + 2(a^ + c^^)]((MB-m5)-(^z,.-^,))^ (5.24) 

Sl""^' = -^-ia"" + a^){{MB-mb)-{En'-E,)f. (5.25) 
2-C/c 

As it was already mentioned, S^^'^ could be obtained from equations (5.20)-(5.24) by 
substituting nic —rric. 



5.3 Real gluons contributions 

The corrections to the hadronic tensor coming from the real gluon emission have also been 
calculated, they could be read off from different works (see for example [f 7], [26], [32], [20]- 
[22], [27]). The corresponding expressions are however pretty lengthy and we are not going 
to use them in our calculation of the corrections to the sum rules. Instead we are going to 
derive approximate formulas for the moments of the structure functions, utilizing the fact 
that the energy release in the decay is of the order of nib — m-c >> Mb — Af^ual ~ ^QCD, 
where M^ual is the hadron invariant mass scale at which the local duality starts. In fact, 
Mb — Afdual is the size of the resonance domain, in which we are going to expand (speaking 
more presicely, we are going to expand in (Mb — Mdual)/ Ec)- 

Let us first define the quantity we are going to calculate. The contribution of real 
gluons to the hadronic tensor I^™* in the leading order in l/nib is defined by the following 
expression, containing c quark and a gluon g in the intermediate state: 
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WJ^' = ^(2vrr^4(pQ-g-px)— (6|4(0)|X>(X|j;(0)|5> 

X = cg 

where A; is the gluon momentum and the sum runs over gluon polarizations. The matrix 
elements (c^|jj,(0)|6) are given by the diagrams on Fig. 5.1. 




Figure 5.1: Feynman diagrams for the real gluon contributions to the hadronic tensor. 
Functions Wj(e, q^) satisfy the following sum rules: 

Z-K JMn-En* J(Mn-mt.)-(En*-Er) 



Let us denote 



' 2tT J(MB-mi)-{EB*-Ec)+\ ^ 

are the contributions of virtual gluons, and 

Tfc,real 1 f j k real 

/, = — / de e" w'^"^', 

ZtT J{MB-mt)-{EB*-Ec)+X 



(5.27) 



^ t«f ^°^(e, q^)e'de. = if + /f + /f (5.28) 

27r J(MB-m5)-(i;B.--Be) 

where are the leading order in moments, 
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are the contributions of the real gluons, Ec = + q^, Cdual = y^dual ~^ ^ ~ y -^D* "I" ^ 

is the excitation energy where the continuum spectrum of hadrons starts. The lower limits 
of integrations correspond to ^ go ^ — Ec — X, with A being the gluon mass needed 
for regularization of the infrared divergency. 

We now need to calculate the right hand sides of Eq. (5.27). To be able to perform 
the calculation we will make one more approximation. In fact in Eq. (5.27) the integration 
goes over the resonance domain of the final hadron energies. We will assume that this 
resonance domain is narrow compared to the maximal energy release in the process. Then 
Cdual can be considered as a small parameter along with all gluon energies in the integration 
domain and we can expand in it. (Note, that in fact we need to calculate moments of 
different structure functions, and not the functions themselves. This allows us to change 
the order of integration and first expand and integrate in e, which makes the calculation 
easier to do). Let us denote 



Let us introduce for convenience the variable v' = —q/Ec, which is the spatial velocity of 
the final quark. Then the results of the calculations look as follows. 
For VV currents: 




yv,o 



{edual [ Ec (4 nibEc-i mbnic - 2Ec^) + 



STTnib E^ 







Rl 



yv,2 





dirE, 



'C 
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(2 m,2i?3 - m, (3 m.W^El + 2 (m,^ + v'^i?,^)))^— ln(i 
2-6^c^(-2 + ^ln(i±^))ln(^)}, 



^2 = ^ ^<^-i[2^-H^1^(T^)]' 



2 97r K 



1 , .1 



(_^2 ^3 ^ (2 ^r'^Et + (3 m,^ + 4 v'^i?^) )) — — ln(- 



2m,^2(2_ in(_^))ln(^)}, 



For FA currents: 



RV^^ = ^{-^[2^.(^.(^.2, 2v'2^2^_^.3^, 



2m6 



i^r^'^ = ^i^[2^.-^ln(i±M) ], 
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For AA currents expressions are obtained from the ones for VV currents by replacing 
nic —Trie. 

In these calculations we did not retain terms of the order (as/'iT)x (nonperturbative 
corrections) to the corresponding moments, because these terms have extra cis/vr and are 
expected to be suppressed in respect to the nonperturbative corrections. In other words, 
here we are not considering perturbative corrections to nonperturbative corrections. 

Note, that first and second moments of different structure functions are simply re- 
lated to each other, which again is a refiection of the heavy quark symmetry: 

(Ec - mc)Rl = 2mbRl = -R\ = RI, (5.42) 

and 

(Ec - mc)Rl = 2mbRl = "^5 = ^l- (5-43) 

Terms violating the heavy quark symmetry show up only in zeroth moments. 

In this section we used notation slightly different from [41] and from the previous sec- 



tion, but the relation between the variables is clear: w = Ec/nic = I/Vl — v'^, y^w"^ — 1 = 
|v'|/-\/l — It is easy to see that the calculated expressions are not singular as v' 0. 
An explicit check shows that the infrared divergency is cancelled. To do this it is 

sufficient to see that in the virtual gluons part the only coefficient that has infrared divergent 

i ■ V A 

terms is a ' : 



In— ... = -- — [^===ln(w-F Vw^ - 1) - l]ln— ... 

TT 4 A O TT \/ w — I ^ 

A: a., rrir. , w ^ \/ w'^ — I rrir. 

= , I In , = - I In— ^ ... 

3 TT ^2\v'\E^ yj-^/w^^ J A 

2 a, r 1,1 + |v'| , , rric 

where we only kept logarithmically divergent terms. 
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Chapter 6 

The QCD constrained model of 

semileptonic decays of the B 
mesons in the heavy quark hmit 

6.1 Introduction. Description of the model 

In this chapter we introduce a model of semileptonic decays of B meson into charmed final 
states B IvXc, based on QCD. In this model we will use experimental information and 
heavy quark symmetry along with the constituent quark model spectroscopy in order to 
identify possible final states and to parametrize the structure functions of the decay. We 
use these states to saturate the phenomenological part of the OPE sum rules (4.9) and then 
use the sum rules to fix the hadronic structure functions. 

Let us start with the description of what is known about the final hadronic states of 
the decay. The lowest energy final hadronic states seen in experiments with B mesons are 
(all D mesons electric charges are zeros, like in decays B~ X^W) : _D(1870), -D*(2007), 
-Di(2420) and _D2(2460) [45]. In this chapter we use the notion of the heavy quark spin 
symmetry [29], [30] and consider final states as doublets of that symmetry. Therefore in 
our model we will consider B meson decaying into leptons and hadronic resonances: pseu- 
doscalar D and vector D* , axial- vector and tensor Di, scalar and axial- vector D**, D^* 
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and pseudoscalar and vector J}***^ D^**. In the constituent quark models the last doublet 
represents the first radial excitation and for some quark- antiqu ark potentials is degenerate 
with the third one. Anything lying above will be considered continuum states. We consider 
the members of the heavy quark doublets to have degenerate masses. We hope that this 
picture provides a good approximation for the final states of the decay and that these final 
states will saturate the OPE sum rules. In the present work we will limit the discrete states 
by the final states listed above, however one can of course add more final states such as the 
first radial excitations, etc. We note, however that the structure functions for the case of 
three heavy doublets in the final state are completely fixed by the first three OPE sum rules 
and therefore could be considered most suitable as the model with the next to leading order 
accuracy. The results obtained in the present section will be used as zero approximation 
for the model including next to leading order corrections to the structure functions (see 
Section 7). 

The resonances listed above could be viewed in the following way in terms of the 
constituent quark model [29], [30]. Since the c quark is heavy, in the first approximation 
we can neglect the interaction of its spin with the light degrees of freedom. Then we can 
classify the states of the light subsystem as 1^5*1/27 1^-^1/27 ^^^3/2 '^^^1/2^ where the 
usual spectroscopic notation n^^'^^lj was employed: n is the radial quantum number, s is 
the spin of the light quark, / is the orbital momentum of the light quark and j is the total 
momentum of light degrees of freedom. The last state is the first radial excitation. Now we 
add the spin of the c quark. Each of the listed states gets doubled thus forming the heavy 
quark doublets and we get the following classification [30] (now in the n^'^~^^Lj notation, S 
is the combined spin of the two quarks inside the meson, L is the orbital momentum of the 
light quark, and / is the spin of the meson or total angular momentum of the bq system): 

pseudoscalar D 11^5*0) and vector D* \l^Si); 



The notation for the last two doublets is in no way standard, we use it here in order to 
merely have notation for these particles. 



pseudoscalar D 



axial- vector Di 



scalar D 
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Since, as it was already mentioned, the last two doublets are usually degenerate in 
quark potential models, we will consider them as one and the same doublet number three. 
Therefore, contribution of the last two doublets in the heavy quark symmetry limit will 
always look like a contribution of a combined doublet. 

As a consequence of the heavy quark mass being large one also can expect that the 
mass splitting inside heavy quark doublets (which is in fact the hyperfine splitting in the 
quark model of hadrons) is smaller then that between different heavy quark doublets (the 
latter is due to the LS coupling in this model). 

In the limit of the exact heavy quark symmetry the hadronic tensor is described by 
just one structure function w(go, 'f') and has the following form (same as for the free quark 
decay) 

W^i, = 47rw(go, g^) (-(Mb - go) gf,u + '^Mb v^v^, + ie^i.c.p'v" q'^ - (v^q^, + v^^q^)), (6.1) 

where we stick to the convention e^i^s _ —£^-^23 = 1. Note that W4 (coefficient in front of 
the structure q^q,^ ) is zero in this approximation. Therefore in this limit it is sufficient to 
write the sum rules for the structure function wi 

wi = 4 TT (Mb - go) w(go, q^)- 

Let us note that for a free quark decay w(qo, q^) = I • 6(M^ — 2MBqo + ^o ~ ~ M^). 

As it was described above, in the sum rules we fix and calculate the following 
moments for the V — A current: 

ilM') = 7r / e'K(go,q%rfe, (6.2) 
27r Jq 

where subscript denotes the leading order in l/rriQ expansion and 

e = Mb- qo- Ed- = Ex - Ed- = + q2 - ^M^. +q2, (6.3) 

emax = Mb-E*d. (6.4) 

The variable e is just the energy of the final hadron counted from the energy of the D* . 
According to the sum rules, these moments of the phenomenological structure functions 
should be equated to the corresponding moments of the theoretically calculated function 
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y^theor^ The number of moments which we are able to calculate is limited by the accuracy 
of the operator product expansion (see discussion in [10] and in Section 4.1). 

If one tries to approximate the final hadronic state by a sum over heavy quark 
doublets and neglects the contribution of continuum states, then for w(go,q^) one can 
write: 

(go, q') = 4 vr ^ (Mb - qo)wHq^) 8{MI - 2MBqo + ql - q' - Mf) 

i 

= 2TrJ2wH^^^)HMB-E,-qo), (6.5) 



where Mi is mass and Ei = \J Mf + is energy of i's heavy quark doublet and the sum 
runs over the heavy doublets. Note that the second ^-function S(MB+Ei — qo) was dropped, 
since it does not contribute to the sum rules, as it was explained in Section 4.1. 

Let us say that both nonperturbative and radiative corrections violate the heavy 
quark symmetry, for them the equation (6.1) does not hold. We will not consider them in 
this section but will take them into account later as small perturbations of the sum rules. 

The function w satisfies the following leading order sum rules: 

i 

J2(E,-En*)w^ = A(l-^), (6.7) 

i 

Y^iE, - En^ f = AHl- ^ f + Y (1 - |f ) (6.8) 

i c c 

As we can see, these equations express phenomenological structure functions through 
quark masses, A and fi'^. If we could solve these equations for w's and express w's through 
A and ^u^, then we could fit A and fi'^ by comparing decay distributions (2.11) with data. 

For the differential distribution in variables Eg, q^, qo we have for the case at hand: 

"^'^ 'V,b\'^w {qo-E,)(2E,MB-q'). (6.9) 



dEgdq^dqQ 2 vr-^ 

Now integrating over go we obtain the expression for the double distribution in q^ 
and Eg-. 

^ = |V.P^(2£. Ms - E ^(Mb - ^. - <e-l«) 
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In the last equation is not independent: 



^2 _ „2 „2 _ M%+Mf + cf- 2MbM, - 2MBqo - 2M,qo 

AMI 



= -q'= " ' ^ — (6.11) 



and 



Ml + M? -q^ , , 

E, = ^ —. (6.12) 

2Mb 



For each resonance the variable q changes between the following limits: 

o 2EJMI - 2MBEe - Mf ) 

» < " < ms-Ie: ■ ("^" 

The distribution in electron energy is now obtained by integrating the double distribution 
over q^: 

2Ee(Ml-2MsEe-Ml) 
ClL / M„-2Ee CI i 



dEe Jo dEe dq^ 

For each resonance in the final state the electron energy lies within the limits 



Mb _ Ml 
2 ^ Ml 



<i?e < ^(1-T7t)- (6-14) 



The general case including perturbative and nonperturbative corrections will be 
considered in Section 7. 

We are now going to consider two simple models for the function w(go,q^), namely 
models with two and three heavy hadronic doublets in the final state of the decay. For these 
models equations (6.6)-(6.8) could be solved and solutions could be used for measuring ji^. 

6.2 Two doublets model 

Let us now try to satisfy these sum rules in the following simple way: assume that there 
exist only two heavy quark doublets with masses Mb,* and Mb,* ■ 
Let us denote for convenience: 

A = Md* - Md*, (6.15) 

= Mjj. - Mjj, . (6.16) 

Note, that in this chapter we do not distinguish between masses of particles within the same 
heavy quark doublet. 
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Now the sum rules (6.6)-(6.8) take the form: 

+w'^ =1 (6.17) 

n? A2 = A(l - + Ed*), (6.18) 

n? iK-f = iK\l-^f + ^{l-^)){Eo*^+Eo* f. (6.19) 

We see that the system of equations is overdetermined and there is no way to satisfy these 
sum rules at arbitrary q^. 

However let us expand the RHS in 



where 



El = ml + q2. (6.21) 



Then for small velocities of the particles in the final states we can cast the sum rules into 
the form: 

+w'^ =1 (6.22) 
a2 = A-^(Mc. +Mc.), (6.23) 

A^ = ^ -^{Md* + Md*)\ (6.24) 

O I-J j-y^ 

The last three equations can be solved: 

A = ^AA, (6.25) 



.,'=1--^. ,6.27) 

The Isgur-Wise p parameter is defined as the slope of the universal Isgur-Wise formfac- 
tor [30], which is related to the slope of w^(q^): v)^ = 1 — (p — |)(v')^, where v' = 



q/-^M|), + q2. From the equation (6.27) we see that p^ is related to A and A in a very 
simple way: 

9 A 1 , , 
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Finally, 



4 ' Ml. + q2 



=l-{p'--) .rt"..,- (6.29) 



Now, with A = 0.55GeV, Md* = 2.009GeV, Md* = 2.462GeV we obtain: 

111 = ^A(Mn. -McO - 0.37GeV, (6.31) 

which is in agreement with the other estimates and experimental data. 

Summarizing, we can say that the first (Bjorken) sum rule along with the second 
(Voloshin) and the third (Bigi-Uraltsev-Vainshtein-Shifman ) sum rules relate A and the 
mass difference A with nl and the hadronic structure function w . Now using Eq. (6.10)- 
(6.14) one can fit nl in dT/dEg from experimental data on the decay. 



6.3 Three doublets model 

It was shown in the previous secsion that in the case of two doublets in the final state the 
sum rules equations could not be solved exactly at any q^. In fact, we did not have enough 
degrees of freedom to satisfy the sum rules. However if one adds the third doublet to the 
model then the sum rules could be solved exactly. In this section we solve the equations 
for the structure functions in the case of three heavy quark doublets in the final state. In 
fact, if we consider masses of the third and fourth doublets degenerate (which happens in 
potential quark models, as it was mentioned in Section 6.1), the structure function w"^ could 
be considered as a sum of contributions of those two doublets. In this sense the proposed 
model in the heavy quark limit takes care of the fourth doblet. But since it is described by 
just tree independend functions, we will still call it three doublets model. 

Let us first introduce the notation. We will call the three doublets di = {D,D*}, 
(^2 = {D2,Di}, and d^ = {D**,Dl*} and use the following definitions: 

A2 = E,,-E,, = ^Ml+q^-^Ml+q\ 

As = E,, - E,, = ^Ml + q2 - ^Ml + q^. (6.33) 
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In this model in the heavy quark limit we use only three doublets, bearing in mind 
that there is a fourth one, which is degenerate with the third, then we do not separate their 
contributions and use generic notation for both of them. 

Now the sum rules for the three particles take the following form. 

+ w"^ + = 1, 

w^A2 + w^A3 = A(l-^), 

w'Al + w'Al = A^(l-if)'+Y(l-|f)- (6.34) 
The last three equations could be solved: 

The three functions w^, w^, w"^ must be positive, which leads to the following in- 
equality: 

A2A < (1 - ^) + ^(1 + ^) < A3A. (6.36) 

Of course, the last inequality is specific for the three resonance doublets case of our 
model, but it still gives some physical limitations on the quantities involved. Let us note 
also that in the limit ^ it looks as follows: 

2 

AmaA < < AmgA, 

which, in turn, could be viewed as limitations on the values of ji^, Am^ and A. For example, 
if one takes Am2 — 0.6GeV, A ~ O.SGeV, then one gets: 

0.45GeV2 < 111, 
0.6GeV < Ams, 

which is in accordance with the estimates for of the works [10] ,[9], [48], [53] (in these 
works the inequality < fj,Q was proven; while fj,Q ~ 0.36GeV^). In fact, however, these 
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inequalities do not take into account perturbative corrections to the sum rules, which could 
somewhat alter them. 

Expanding at small we find 



2A3A2 3A3A2^' ^ ^ 

where = q^/(m^ + q^). From the last expression we can see that the Isgur-Wise p- 
parameter in this case is 

2 I ^ ml A(Am3 + Am2) pi 
^ 4 Ml, 2Am3Am2 3Am3Am2'^' ^ ' ' 

where Am2 = M^^ — and Am3 = M^^g — M^^ and we have taken into account the fact 
that the /^-parameter is defined for hadrons and not for quarks. 

For the differential distributions one can again use formulas (6.9)-(6.14). 
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Chapter 7 



Model of the decay including 
perturbative and nonperturbative 
corrections to the sum rules 

7.1 Description of the model. 

In the previous chapters we have considered the model of semileptonic decays of heavy 
flavors in the heavy quark limit. We considered the final states resonances to be heavy 
quark doublets and also neglected terms breaking the heavy quark symmetry in the sum 
rules. In this chapter we are going to take into account the fact that the resonance doublets 
masses are split and that the hadronic tensor and, therefore, the sum rules, contain terms 
violating the heavy quark symmetry. We are going to use the previous chapters results as 
a zero approximation model and consider perturbative and nonperturbative corrections as 
small perturbations to it. 

Let us consider a model with n heavy quark doublets. To accommodate deviations 
from the heavy quark symmetry we now need to consider mass splitting within the heavy 
quark doublets as well as independence of different structure functions in the hadronic 
tensor (they were not independent in the heavy quark limit, see Eq. (6.1)). 

Now we have 2n resonances with different masses, we will consider mass splittings 
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of the doublets to be of order I/ttiq. 

For the masses of the first four resonances experimental data is available, and there- 
fore we do not need to consider them as parameters to be fixed from the sum rules [45]. 
For higher resonances, however, no data exists, but as numerical analysis of the model will 
show, they are not going to play a major role for their contribution is going to be small 
anyway. 

The hadronic tensor now has all structure functions independent from each other. 
Three of these functions contribute to semileptonic decays. The operator product expansion 
was calculated only with the accuracy I/itiq. This limits the order of the sum rules we can 
use to restrict the structure functions by three. It is valid to include l/rriQ corrections to 
the first two sum rules while in order to calculate them for the third one we would have to 
know the higher order terms of the OPE. 

First let us introduce the notation. Since now we are considering deviations from 
the heavy quark symmetry, we assign each resonance its own distinctive mass Mi, where 
for odd i's Mi, Mj_|_i come from the same heavy quark doublet. Let us also denote Wj, Wj"*"^ 
the corresponding structure functions, where j stands for the structure function number 
and i,i + 1 denote the functions coming from the same doublet for odd i's. Note that 
functions refer to the separate resonances within the heavy quark doublets, while the 
previously defined functions which did not have subscripts referred to the contributions 
of the whole doublets. Let us also denote the leading order moments as I^j. In this section 
we only consider V — A current. 

Our model is now formulated as follows: 



where Mcont is the hadronic invariant mass at which the continuum spectrum starts. 

In our model we will assume that Mcont = Afduab other words, we assume that local 
duality starts at the same scale as continuum spectrum. Then in the LHS of the sum rules 
only resonance contributions should be counted. 



phen/ 2\ 

% (go,q ) 



27r ^ «;;.(q2) S{Mb - .Jm^ + q2 - go) + 
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The first three sum rules now read: 

E K + <') = ^, (7.2) 

8 = 1,3,5,... 

J2 w){E,-ED') + wf\E,+r-ED') = i], (7.3) 

8 = 1,3,5,... 

E w](E, - ED*f + wf\E,+i - ED*f = if. (7.4) 

8 = 1,3,5,... 

In the last equations Ij are the moments including all perturbative and nonperturbative 
corrections found in previous sections. 

We are interested in the lowest order in l/rriQ corrections to the sum rules. Let us 
expand energies of the final particles in heavy mass splittings: 

E. « - %^ - - if- - - 

Then we have for the sum rules: 

E (w) + <^)(^Hi - En^ f - ^;.^%i^(^.+r - En^ f-' = /^ (7.6) 

8 = 1,3,5,... ^' + 1 

For w^j + '^]'^^ we can write 

w] + w'+^ = (w] + w'+^)o + 6w] + 6w'+\ (7.7) 

where (wj)o are the leading order structure functions which are determined from the sum 
rules in the heavy quark limit (6.6)-(6.8), and Sw'j are corresponding corrections. In fact, 
by comparing (7.1) and (6.5) one finds: 

= {w\ + w\+^)o, i = 1,3,5. 

Let us take into account the fact that SMi ~ ^qcd/'^Q therefore Wj"*"^ in the 
second term of Eq. (7.6) could be taken in the leading approximation. By the virture of the 
heavy quark symmetry both (wj)o and [w^-'^^)o could be expressed through the universal 
Isgur-Wise formfactors [30]. The explicit expressions for these functions are given in the 
appendix C. 
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Since it is only the sum of corrections for the members of a heavy quark doublet 
which enters the sum rules and the expressions for lepton energy distributions, we will 
introduce a special notation for it: 



where Iq j are defined in Eq. (6.2). 

This system of equations relates corrections to the structure functions to the mass 
splittings inside the heavy doublets. For the model to be self consistent in the zeroth 
approximation we take resonance doublets masses equal to the masses of heavier particles 
of the doublets. This does not matter in the zeroth approximation itself (heavy quark 
symmetry limit) but does matter when we consider nonperturbative corrections to the sum 
rules coming from the mass splittings within the doublets. Again, for the first two doublets 
the mass splittings are known experimentally [45]. 

We see that the OPE sum rules provide the system of equations for the functions 
^Wj, which is determined to the same degree as the system of equations for the functions 
(wj)o. In fact for V — A current we have three systems of equations, independent for each 
structure function number. 

7.2 Three doublets model 

For three heavy doublets in the final state the system of equations (7.8) can be easily solved, 
for it has the following generic form: 



SW; = Sw) + Sv/+\ i = 1,3,5. 



Finally corrections to the sum rules can be written in the form: 



8 = 1,3,5,... 




X + y + z = a. 



yAl + zAl 



c. 



The solution is: 



X 



c - 6(A3 + As) + aAgAs 
A3A5 
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_ 5A5 - c 
y - AsAs-Ai' 

&A3 - c 
' - AsAs-Af 

Let us denote: 



Ai = -Ei - -Ed*, 

^/^. = - /q^^. - ^ {wr^ )o ^. A 



8 = 1,3,5, 



Then for each structure function the equations are: 



and the solution is given by: 



bW] 



W/As + W/As = bl], 

bwf^i + bw|^l = bi], 



bif - blj{As + As) + ^/°A3A5 



AaAs 



bW- 



bijA, - bij 



' A3A5 - A2 ' 

bnA^ - bif 

7.3 Differential distributions for the model including pertur- 
bative and nonperturbative corrections 

One of possible uses of the proposed model is to describe the lepton energy and invariant 
mass distribution in the semileptonic decays of heavy mesons. In this section we will show 
that the information about structure functions derived from the corrected sum rules is 
sufficient in order to write expressions for the differential distributions. Unlike the pure 
OPE results, this model predicts smooth distributions, which do not have anything like 
^-functions at the ends, and still is totally based on OPE. 

The hadronic structure functions of the model are given by the following expression 
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wf^\qo, q2) e{Ml - 2MBqo + - MI^) 
27r ^ w]{q^)2E^6{Ml - 2MBqo + - M, 



2^ 



w'^°^\qo, q^) e{Ml - 2MBqo + q' - M^^). (7.10) 

In the last equations we have taken into account the relation 

/ dq^dqofiqo, q^)2E,8{Ml - 2MBqo + q^ - Mf) 

J decay phase space 

= I dq'dqof(qo.q'o-q')6(MB-qo-jMf+^)^ (7-11) 

J decay phase space 

with corresponding limits of integration. Because, it is the double distribution in Ee,q^, 
that we studied in the previous chapters and experimentalists are going to measure, in what 
follows we will work in variables (go, q^)- The form of the equations in the variables (go, q^) 
could be readily restored. 

The full differential distribution is given by the Eq. (2.9), which we rewrite here 
showing explicitly contributions of different states: 

dEedq^dqo ' 32 vr'*'^ 

27r J2 [2g'^i + [4 Eeiqo - E,) - q^] w\ + 2 q\2 E, - go) ^^3] X 

i 

2E,6{Ml - '^MBqo + q^ - Mf) + 

[2q^wl°^' + [4 EMo - Ee) - q'] w^^' + 2 q\2 E, - go) w^'] X 
9{Ml - 2MBqo + q'- M^.J}. (7.12) 

The total width is then given by the formulae: 
d^T 



dEedq dqo- 
$ dEe dq'^dqo 

KhPCl^ v-^ f ,^ , o E, 



,o^i27r^ / dEM^-^[iq^w{ + [AE,{MB - E, - E,) - q^]w\ + 
2q\2Ee-MB+E,)wi,] + 

[ rf<i)[2g2«;™-t + [4 E^iqo - E,) - q^] u;™"* + 2 g2(2 E, - go) u;^*] X 
e{Ml - 2MBqo + q'- M^.J}. (7.13) 
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where <i>,<i>j's are volumes in the phase space in corresponding variables available to the 
particles in the final state, and the sum runs over all resonances. Note that the phase 
spaces available to heavier particles in the final state are smaller then those available to the 
lighter ones, therefore ^i+i > ^i, i = 1,3,5. 

Let us examine the resonance contributions to the total width. First let us denote: 



r = -^[2q^w\ + [4E,(MB-E,-E,)-q^]wl + 2q^(2E,-MB + Ei)wi,], (7.14) 
Mb 

where 

^ ^ M| + - 
2Mb 

For a contribution of a heavy doublet, up to the first order in the heavy doublet mass 
splitting, we have: 



[ dEedq^f+ [ dEedq^f+^ 



dE,dq\r)o+ / dEJqHr+ / rf^eV(r+')o + / dEdqHr+^ 



dEJq'{r)o+ dE,dq'Sr+ dEJq'{r+')o + 

dEJq\f)o+ [ dE,dq^Sf+' 
dEJq\r)o + ir+')o + iSr + Sr+')]+ [ dEJq^iDo, (7.15) 

where (/')o are taken in the leading order (heavy quark symmetry limit) and SP are correc- 
tion to them. From the last equation it follows that there are regions of the decay phase space 
where only one resonance of a doublet is present. Therefore unlike heavy quark symmetry 
case, now we need to know not only leading order structure functions describing resonance 
doublets, but also the leading order contributions of individual resonances. These contri- 
butions could be expressed through the universal Isgur-Wise form factors [29], [30], which is 
done in Appendix C. 

In the three doublets model, all the integrands in the last line of Eq. (7.15) are 
known. To find (/')o + (/''''^)o we can use solutions (6.35) of the model without corrections. 
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and to find (SP + SP~^^) we can use (7.9): 

Sr + Sr+' = -^[2qM + (4 E,(Mb - - E,) - q^) w\ + 2 q\2 - Mb + E,) wi,] - 
Mb 

{^[2q'w{ + (4 E,(Mb - E, - E,) - q') w\ + 2 q\2 E, - Mb + E,) u;^ }o + 
Mb 

^[2q'w\+^ + (4 E,(Mb - ^,+1 - Ee) - q') w'+^ + 2 q\2 E, - Mb + ^3+'] " 

Mb 

{^[2q^w\+^ + (4E,(Mb - ^,+1 - Ee) - + 2q^(2E, - Mb + ^,+i)^+']}o 

Mb 

^'-[2qHw{ + (4E,(Mb - E, - E,) - q^)6w\ + 2q^(2E, - Mb + E,) Sv/^] + 



Mb' 

[2qHw\+^ + (4 E,(Mb - ^,+1 - Ee) - q^) Sw\+^ + 2 qH2 E, - Mb + ^^+'] 



-£^8 + 1 

Mb 



^ ^[2q\Sw\ + Sw\+') + {4E,{Mb - E, - E,) - q^) (bw\ + bw\+^) + 
Mb 

2 q\2 E,-Mb + Ei){ bwl, + bwi+^)] 

= ^[2qHW{ + (4 E,{Mb - E, - E,) - q^) bW\ + 2 q^{2 E,-Mb\ E,) bW^]. (7.16) 
Mb 

Finally, we get for the total width: 

r'"^ r 1 



r = \Vch\''i^A I dEJq^{4{2E,MB 



Mb 

E 



(Mb - E, - E,) 



Mb 



[2q'bW{ + (4 E,(Mb - E, - E,) - q') bW\ + 2 q\2 E, - Mb + E,) bW^]} 



+27r J2 I dE,dq^[2q^iw{)o + iiE,iMB-E,-E,)-q^)iw'^)o 

, = 1,3,5 J'^'-'^' + i 



+2qH2E,-MB+Ei)(w';,)o] 



Mb 

+ [ dEJq^dqo2Tr[2q^w'r' + (4 E,{qo - E,) - q^) u;™"* + 2 q\2 E, - go) w^^*] 
xe{Ml-2MBqo + q^ -Ml^)}. 

In the last equation the first term is the leading order contribution, it contains functions 
defined in Eq. (5.27)) which are found by solving leading order sum rules equations (6.35), 
second and third terms are the nonperturbative corrections and the fourth is perturbative 
corrections to the total width. The functions (wj)o are read off from the equations (C.2), 
(C.ll), (C.20), go = {M% + M| - q^)/2MB, where Mx is the mass of the corresponding 
particle. Numerical analysis shows however, that contributions of the phase space regions 
— $8-1-1 are negligibly small. 
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The energy distribution is the integrand of the integral over Eg, and the double 
distribution in Eg and is given by the integrand of the integral over dEgdq^ . For each 
resonance the variables change in the following limits: 

2 2E,{Ml - 2MbE, - Mf ) 



< < 



Mb - 2E, 
Mb,, MJ 

As for the contribution of continuum states, we recall that in our model it was chosen 



< < —(1 - ttt)- (7-17) 

2 Ml' ^ ' 



to be the same as radiative gluon emission correction to the free quark decay starting at 
some invariant hadronic mass M^^^. These corrections have been calculated many times 
before (see for example [17],[32],[20]-[22], [26], [27]), and we are not going to reproduce them 
here. The range of the variables for the continuum is as follows: 

Ee + < ?o < 



o 2EJml - 2mbEe - Ml- ) 
< g2 < — simZ^ (7.18) 

m - 2Ee 

^ nih , Ml- , 
< < -^(l--jr). (7.19) 

Use of the quark variables is dictated by the nature of the calculation in which one 
considers decay of a quark into a quark and gluon. However use of hadronic variables (which 
is effectively use of Mb instead of m^) is also permissible within our accuracy since we did 
not take into account corrections of the order cisAg^j^/m^. 

This effectively concludes formulation of our model. 



7.4 Numerical analysis. Importance of the nonperturbative 
corrections to the sum rules. 

In this section we are going to investigate how important the nonperturbative corrections 
to the sum rules are from the point view of affecting predictions of our model. This section 
is mostly illustrative. More work is needed in order to make reliable connections between 
the proposed model and real experiments. 

Let us make some assumptions about different quantities that enter the model. First 
about hadron masses. For them we take the following values [45]. For i?-meson mass 
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0.05 0.1 0.15 0.2 



Figure 7.1: Hadronic structure functions w'(q^), i = 1,2,3, for the model without nonper- 
turbative corrections to the sum rules. 




0.05 0.1 0.15 0.2 



Figure 7.2: Hadronic structure functions (w\-\-w^^^)(q^) for the model including nonpertur- 
bative corrections to the sum rules. These should be compared directly with the functions 
on Fig. 7.1 
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Figure 7.3: Functions p(Ee, q^)o of the resonance doublets for the model without nonper- 
turbative corrections to the sum rules for x = 2Ee/MB = 0.2 . 
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Mb = 5.2GeV. For the first resonance doublet: Md = 1.87GeV, Md* = 2.01GeV; for the 
second: M^* = 2.42GeV, M^^ = 2.46GeV. About the third resonance doublet experimen- 
tal information is not available and for testing of our model we take the following values: 
Md** = 3.8GeV, M^** = 3.83GeV. For the matrix elements of the operators in the OPE we 
can use the following information. First, as we mentioned earlier, fj,Q = .36GeV^ is known 
from the hyperfine mass splitting between B,B*. For there is no direct measurements, 
but we can use the inequality fi'^ > /Jq to get an idea about its possible range. For our 
estimates we take fi^ = .4GeV^, which is consistent with the other people's estimates (see, 
for example [23]). 

We did calculations of our model both in the heavy quark symmetry limit (without 
nonperturbative corrections to the sum rules), and without heavy quark symmetry (nonper- 
turbative corrections to the sum rules taken into account). We calculated hadronic structure 
functions, which in fact determine the corresponding hadronic formfactors, functions /' 's 
(7.14), which represent contributions of different resonances to the differential distributions, 
and electron energy distributions. 

The hadronic structure functions w'(q^) and w^(q^) are plotted on Fig. 7.1 and Fig. 
7.2 correspondingly. One can see that the functions are significantly different, which signals 
that nonperturbative corrections to the sumrules are relatively large and important. Same 
conclusion follows from the analysis of plots of functions /' 's on Fig. 7.3 and Fig. 7.4. The 
most different between symmetric and non-symmetric cases is behavior of contributions of 
the first and second resonance doublets, while the contributions of the third one in both 
cases is relatively insignificant. This could be an indication of the fact that two doublets 
approximation of Section 6.2 is sufficient (with nonperturbative and perturbative corrections 
to the sum rules taken into account). 

On Fig. 7.5three normalized electron energy distributions are plotted. Two of them 
- for the cases of the model with and without nonperturbative corrections to the sum rules 
taken into account. The third curve shows the electron energy distribution for the free quark 
decay. We plotted all of them on the same graph in order to demonstrate differences in the 
electron spectrums. We see that the model predicted spectrum is different from the free 
quark one. Let us also mention, that despite the fact that the electron spectrums for the 
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Figure 7.5: Normalized electron spectrum plotted for the model with and without non- 
perturbative corrections. For comparison, electron spectrum for a free quark decay is also 
show. 

models with and without the corrections seem close, the determination of phenomenological 
parameters from experimental data based on the model with and without the corrections 
of course will be different. 

7.5 Summary and outlook 

The proposed model of the semileptonic decays B XJi^ is based on some general heavy 
quark spectroscopy arguments and the operator product expansion sum rules. In this sense 
the model is so to speak "model-independent", for it does not utilize any quark- antiqu ark 
potential nor it introduces any phenomenological parameters like "Fermi momentum" which 
are not coming from the QCD based analysis. All the parameters of the model have a clear 
meaning in QCD, like hadronic matrix elements ^u^, /Jq, A, etc . 

There are still some unanswered questions about the model, however. One of them 
is about the fourth resonance doublet, which in potential quark models is degenerate with 
the third one. As it was shown, this is not a problem for the heavy quark symmetry case, 
but the issue has to be addressed in the full model as well. One of the possible solutions 
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could be considering more sum rules, not just for the V — A current but also for V and A 
currents separately. This would provide more sum rules and could describe the structure 
of the final state in more details. Another approach could be not using the notion of the 
third and fourth resonances at all, but instead use the sum rules to restrain part of the 
continuum spectrum in the domain where local duality does not work (one can think about 
it as of some broad resonance merging with the continuum spectrum). 

The other direction of improving the model and making it possibly more accurate, 
is to use Borel transformed sum rules. This would effectively suppress contributions of 
the higher states along with continuum states, and at the same time better describe the 
contributions of the lowest resonances. 

Numerical analysis of the effects of perturbative corrections also needs to be done. 

The main question about the model has to do with the accuracy of the operator 
product exopansion itself. We have seen that the nonerturbative corrections of order I/ttiq 
to the OPE sum rules strongly affect predictions for the structure functions of the reso- 
nances. What would happen if the next order corrections were taken into account? How 
this would affect determination of the value of fi'^ from experiment? The I/ttiq corrections 
have been considered in recent works [12], [25], but there is still no conclusive numerical 
estimates for hasronic invariant functions, which could be used in the considered model. 



75 



Chapter 8 



Conclusions. 



In the first part of the work a model independent operator product expansion approach was 
used to calculate differential distributions in semileptonic decays of heavy fiavors in QCD. 
Double distributions in lepton energy and invariant mass of the lepton pair for the massless 
lepton (electron) in the final state and the energy spectrum for massive r-lepton have been 
calculated. 

Based on expressions for hadronic invariant functions, the full set of sum rules for 
heavy fiavor transitions up to order I/ttiq was calculated along with perturbative corrections 
to them up to order asAqcD/mQ- 

A new model of semileptonic decays of i?-mesons was proposed, which is based on 
the model-independent OPE sum rules approach. It was shown that the OPE sum rules 
provide sufficient number of restrictions to fully fix contributions of different resonances to 
the hadronic tensor. In fact, in this approach it is possible to calculate different hadronic 
structure functions for the transitions induced by the weak currents. Important here is 
that the proposed model is so to speak "model independent". It does not use any potential 
quark model calculations and is entirely based on very general symmetry and spectroscopic 
considerations. 

The proposed model could have several uses. On one hand it can realistically de- 
scribe the end point (maximal lepton energy) part of the lepton spectrum as well as double 
differential distributions in lepton energy and combined invariant mass of the lepton pair. 
On the other hand, it provides relations between the lepton spectrum and important phe- 
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nomenological parameters describing strong interactions in heavy mesons, which are given 
by the matrix elements of operators in the operator product expansion. This should enable 
experimentalists to extract the OPE matrix elements from the experimentally measured 
lepton spectrum. 
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Appendix A 

The phase space of the decay 



In this Appendix we briefly outline the derivation of the Lorentz Invariant Phase Space 
(LIPS) for the case of the inclusive semileptonic decay Hf, tvX . Let us introduce the 
invariant kinematical variables in the following way. Let P be the 4-momentum of the 
decaying particle, so that Mjj^ = P"^ q = -\- \s the 4-momentum of the lepton pair, 
jjt and are the momentums of the emitted charged lepton and neutrino correspondingly. 
We introduce m\ = p^, - the invariant mass squared of the born hadronic state and 
= Pxfy, ^ the combined mass of the hadrons and z/. We have the following relations: 
px = P — q and pxi? = P — Pr- The introduced invariant variables are related to the ones 
used in [13] in the following way: 

m\= Mjj^ + q^ -2MH,qo. (A.l) 
= M^^ + - 2 Mh, Er. (A.2) 

The inclusive decay rate in the normalization of [13] is given by the following ex- 
pression: 

I X max 9i/9 99\i9/ 

1' = 100 6 / , dm'x\M(m'x^,m'x,q')\U^P,px,Pr,P^), (A.3) 

where | (m^-, m^, is the matrix element which describes the transition into the 
hadronic states with mass it is given by the expressions (2)-(8) in [13]. The lower 

boundary of is given by the mass squared of the lightest flnal hadronic state possible 
in the decay (_D-meson mass for h c and pion mass for h u transitions ), while the 
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upper boundary of is {Mhi, — m^)^. The d^{P,px-,VT-,Vi^) is the three particle LIPS 
for a particle with the 4-momentum P going into three particles with 4-momentums px^Pr 
and p^. It is given by the formulae: 

where 9(a) denotes 9(ao), and uq is the zero component of the 4-vector a. For the three- 
particle phase space one can write: 

/dm?' 
{P,pxi^,Pr) {pxi^,Px,Pi^) Km\^ - p\^) (A.5) 

where d^ (a, 5, c) is a two particle LIPS for a particle with 4-momentum a going into particles 
with 4-momentums h and c. The equation (A.5 gives the decomposition of a 3-particle LIPS 
into 2-particle LIPSes. 

In order to introduce a new variable into the phase space integral one can multiply 
the equation (A.5) by "1": 

'' dq^8{q^ -{p,+p,f )e{qo) = l. (A.6) 

Performing the integrations over all variables except for m\,q^ and m\-, which amounts 
to integrating the delta functions in LIPSes along with determining the conditions for the 
delta functions to be non-zeros, we arrive to the following formulae: 

1 [(^H^-'mr? rm\^ /-SmaxC^X'^xJ 2| / 2 2 2m2 

(A.7) 



where 



1 M"^ — 

2/2 2 \ r i\j2 I 2 2,2 -fffc t 2 / A o^ 

<lrmn(mx,mxi,) = ^ [ + "^X " "^X + "^r "^X (A.8) 

^ "*Xi^ 



(A.9) 



Zi Til ^ 

^A(M^^,m2,m^-) (m^- - m^) 



A.ll 
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and 

A(a, b,c) = + + -2ab-2ac-2bc. (A. 12) 

The physical meaning of the function A(a, 5, c) lies in its relation to the square of the spatial 
momentum p*^ of particles born in a two - body decay in the center of mass reference 
frame, for example: 

*2 _ *2 _ 

The fully differential distribution (3.1) in the invariant variables looks as follows: 

^ \Mim\,m\,,dq'f. (A.14) 



Px,-Pr - • (A.13) 
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Appendix B 



Hadronic invariant functions 



Here we present the results of calculations of different hadronic invariant functions hi, 
introduced by eqs.(2.15), (2.42). The structure functions Wi are simply related to hi by 
eqs.(2.16) and (2.52). We use the following notation: qq = q ■ v, ^ = Qq — and z = 
ml-2mQqo + q^-ml. 

For the Vector X Vector functions we have: 

/^r = -[(-Q--,-go)-(M^-M^)^(^+^)];- 

— [ ^ ( (4 - 3go)( mq - mg - qo) + 2q^) + 

A (go {mq - - go) - ^ ) ] ^ - ^ {mq - nig - qo) ^ , (B.l) 



2 18 1 

-[2fil;{mq - mg) - 5 fil; qo + 7 fil qo] ^ - -mqfilq^—, (B.2) 

hr = 0, (B.3) 

hV = -^{j^l-l^%)^, (B.4) 



3mQ 

hr = l-li^^(^o-^l)-4^l)]^ + I^U'j,. (B.5) 

To get the functions hf^ for Axial X Axial tensor from hj^ one should substitute nig 
by (— nig) in eqs.(B.l - B.5). 
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For the Axial X Vector tensor only one invariant structure survives: 

l,VA l|ro2|5.2 2\?0-|l ,^22^ 



Summing up we get the result for the full hadronic tensor h 
= - [ 2 (mg - go) - {p% - ^l) ] ^- 

' (4 + 2 - 7 mg go + 3 go) + (4 go {mg - go) 



10 , 2 2 M 1 r28 2 2 /8 20 

- [ 4 mg + (//^ - Mg) J - - [y M,r go + ( 3 "^Q - y 

16 2 2 1 



1 



r/i2|10/2 2\"iU-iJ- "2 2 

[4mg + y (m^ -Mg) — ] " 



go T 1_ _ 8 2 2 J_ 

2 q q ^3 



/l4 



3mQ 



mq' 3 
(m^ - Mg) 3 7 



1 2 
I ~ 3 



2x ?0 



1 8 



[5(M^-M^)f^-4//^)]- + -//2 q2 



y3 • 



Appendix C 

Structure functions of the 
resonances in the heavy quark hmit 

In this appendix we are going to write down expressions for the contributions of separate 
resonances to the structure functions in the leading order in the heavy quark mass. We 
need these expressions for the formulation of the full model of the decays which takes into 
account independence of different structure functions and mass differences between particles 
in the same heavy quark doublets. Due to these mass differences there are regions in the 
phase space of the decay, where only one (the heavier) resonance of the doublet contributes. 

These formulas could be obtained using the results of the original works of Isgur and 
Wise [29], [30]. 

Now we need to write down parts of the hadronic tensor coming from each resonance 
in the heavy quark limit. For our purposes we only need to know the structure functions 
for the V — A currents. 

W,, = J2i'^7r)U\pB-q-px)^{B\jlmPx,^){Px,MMm) 

X 

= 27rJ2 KMb -qo- ^Mf + q^) ^ ^{(^I^JWI^ A>(^, q, A|K(0)|i?> 

1 A 

-{B\v}{0)\i,q,\){i,q,\\A,{0)\B)-{B\Al{0)\z,q,\){z,q,\\V,{0)\B) 
+ (i?l4(0)|^,q,A)(^,q,A|A,(0)|i?)} 
= 2TrJ2 w^UmMB - go - + q2), 



where i runs over all resonances of our model (not just the resonance doublets), A runs over 
all polarizations of the corresponding particle and 

A 

-(i?l4(0)K, q, A)(^, q, A|K(0)|i?> + (i?|4(0)K, q, A>(^, q, A|A,(0)|i?)}. 
Now we define the structure functions of the resonances in the following way: 

WfMuii) = -w{g^^ + vj2Vi_,v^ -i e^^^afS ■v"q'^ + w\ + (q^ + q^ v^), 
so that 

w,(qo,q^) = 27r J2w)(q^)6(MB - ^/m~T^ - qo). 

i 

To calculate the hadronic tensor we will also need formulas for sums over spins of 
the vector and tensor particles in the final states: 

^e^(w', A)ei,(w', A) = -g^^ + v'^vl, 
A 

A 16 

where the projection operator P^j, = ^^j, — v'^v'^. 

Now we can write down the contribution of each resonance to the hadronic tensor 
in the heavy quark limit (let us remind the spectroscopic notation: ra^'^+^ij, here n is the 
radial quantum number, S is the combined spin of the two quarks inside the meson, L is 
the orbital momentum of the light quark and / is the spin of the meson): 
for pseudoscalar D 11^5*0) and vector D* \l^Si): 



{D\V^\B) = y^r]^iw)iv^ + v'^), 

{D\A^\B) = 0, 

{D*\A^\B) = y^2M^r]iiw){iw+l)e^-e-vv'^}, 

{D*\V^\B) = i./2M^r,r{w)e^Sap€'v"v'f, 

Wt^AD) = r]f{w){vi,v^ + (vi^vl + v'^v^) + v'^vl}, 

= vIM^UMb + Mx fv^v^ - {Mb + Mx){q^v^ + q^v^) + g^gj, 
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,{D*) = i]l{w){-2w{l + w)g^^-v^v^-2i{l + w)e^^apv"v'^ 

= v!H^{-'^(Mb - qo){MB + Mx- qo)9^, 

+ (3M'^ + 2MbMx - Mx^ - 4MBqo)v^v, + 2i(MB + Mx - qo^^Sapv^qf^ 
-{Mb + Mx - 2qo)(q^Vi_, + q^v^,) - qf,q^}, (C.l) 



2(Mb + Mx- go) 2. ,s ^ 

+2MBVf,v^ + i€f,^apv"q^ - [q^v,^ + q^v^)}, (C.2) 



w = 2{w + l)r]^{w) = — r]^{w) 

wl = 0, 

^2 = Vf{w] 



{MB + Mxf 
M\ 



wl = 0, 



.^2( AMB-qo){MB + Mx-qo) 
,^2 ■ 



Mi 



W2 = Viiw 



i^Ml + 2MbMx -Ml- AMBqo) 
Ml 



^^2. .{Mb±Mx^ 



for tensor \1^P2) and axial-vector Di vf + Jh^^Pi] 



{D;\V^\B) = tV2M^V3v2{w)e^apse"''vyv", 
{D;\A^\B) = V2M^V3v2(w){(w + l)e^^v" -v'f^pv^vf"}, 



{D,\V,\B) = y2MB— ,^2H{-(^'-l)e^ + (-3^;^ + (t«-2)^;;)e.^;}, 
{DM^\B) = -V2M^^V2(w)(w+l)e^^pse"vfv'\ 



C.3) 
C.4) 
C.5) 
C.6) 
C.7) 

C.8) 

C.9) 



+ (1 + w)iw - 2)v^v^ - i{w - 1)(1 + w fe^^c.pv"v'^ 
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+ (1 + w){-2 + 2w + w^){v^vl + vlv^) + (1 + w){w - 2)vlvl} 
= vlH ^^''^^f~'^'h -(MB -Mx- qo){MB + Mx- qo){MB - qo)9, 
+(3M| + 2M|Mx 

-2MbMI - 2Mi - 5M|go - 4MBMxqo - Mj-go + 2MBqo)v^v, 
+i(MB +Mx- qo)(MB -Mx- go)W/3^"/ 
-(2M| - 2Ml - 3MBqo - 2Mxqo + ql){q.v^ + q^v,) 
+ (Mb - 2Mx - qo)q^q.}. 



-(1 + w){w - 2)v^v, - 3i{w - 1)(1 + w fe^,apv"v'^ 
+(1 + w){-2 -2w + 3w^){v,vl + vlv^) -{w- 2)(1 + 



Vliw y^'^^f ^°^ {-3(Mg -Mx- qo)(MB + Mx - qo)(MB - qo)9,. 
+(5M| - 2M|Mx 



-bMBMl + 2Ml - llM|go + ^MBMxqo + Mj-go + QMBqo)v^v, 
+3i{MB + Mx- qo){MB - Mx - go)W/3^"/ 
-(2M| - 2Ml - SMBqo + 2Mxqo + ^oKl^v^ + q^v,) 

-{Mb - 2Mx - go)g^g.}, (C.lO) 

.(i?;) + w^^iDi) = A{w - l){w + lfri{w){-wg^^ - ie^^afjv"v'^ + {v^v'^^ + v^v'^)} 

4{Mb + Mx- qo)\MB -Mx- go) 



M'x 



+2MbVi_,v^ + ie^i,c,f3'v"q^ - {q^v^ + g^Wj,)}, (CH) 



A{Mb + Mx- qofjMB -Mx- go) 
Ml 



= 4{w-l){w+iyr^i{w)= V - — ^y^^v ^^^^^(^). (C.12) 



vlM ^^^'^fj '^'^\ mb -Mx- qoKMB - go), (C.13) 

vjI = rili^)i^^±M^^i2Ml + 2MlMx-2MBMl-2Ml- 

bMlqa - 4MBMxqo - Mj-go + 2MBg2), (C.14) 
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^3 = (MB-Mx-qo), (C.15) 

< = ^riliyj) ^^'' ^ fff " -Mx- QoKMb - go), (C.16) 

AMeMxqo + Mjrqo + GMeq^o), (C.17) 



W3 = 3?72(w) (Mb - Mx - go)- (C.18) 



for scalar D** \l^Po) and axial- vector D^* - 

{D**\V^\B) = 0, 

= -2y2MB7?3HK - <), 

(i^ri^^li?) = 2y2MB7?3H{(t«-l)e^-e.^;^;;}, 

{Dl*\A^\B) = 2iy^2M^r]3iw)e^SafjesVfjv's, 

w^AD**) = ^Uw){-{vy, + vlv,) + v,v, + v'^vl} 

= 4??|(w)^{(Mb - Mxfvi^v^ - (Mb - Mx){q^v^ + q^v^ + g^gj, 

w^i.(i^r) = 4?7|(w){-2w(w - l)£r^i, - v^v^ - 2i(w - l)e^,„ai3v"v''^ 



4%'(^)i;7r{-2(^B -Mx- qoKMB - qo)9, 



Mx 



+(-(Mb + MxY + AMb{Mb - go)) 
+2i(MB -Mx- go)W/3^"/ 

-(Mb - Mx - 2qo)(q^Vi_, + q^v^,) - qf,q^}, (C.19) 



'-V3{'^){-{Mb - go)6'p 



8(Mb -Mx- go) 2, 



+2MBVf,v^ + i€f,^apv"qf^ - [q^v,^ + q^v^}, (C.20) 



3 n 2( ^ 8(Mb -Mx- go) 2( ^ 

wl = 0, (C.21) 



87 



wl = 0, (C.23) 



2, ^^ {Mb - Mx - qo ){MB - go) 
Ml 



8viM^ (C-24) 



2, , -(Mb + Mx)' + 4 Mb{Mb - go) 
Ml 



< = 4y?3» ^ " ' (C.25) 



2. , Mg-Mx-go 
W3 = 8t]s(w) . (C.26) 



For the fourth doublet which is (pseudoscalar, vector), like the first one, the formu- 
las are exactly the same as corresponding expressions for the first one, with the evident 
substitution i]i ^ r]4. 

In the above equations r]i(w), i]2(w) and i]3(w) are the universal Isur-Wise form- 
factors for the three resonance doublets respectively; Mx are the masses of the corre- 
sponding resonance doublets; v,v' are velocities of the initial and final heavy mesons 
(v' = (Mb - go)/Mx), w = v ■ v' = Ex/Mx = (Mb - go)/Mx, (note that w > 1); 
Cj,, e„i/ are the polarization vectors and tensors for the corresponding (pseudo)vector and 
tensor mesons in the final state. 

Functions , w^, w"^ are known from the solutions (6.35) to the sum rules equations 
(6.6-6.8), which in turn determines the form factors r]i,r]2,i]3. The structure functions for 
each of the resonances could be read off from the above equations. 
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